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Abstract 

We study the gravitational Dirichlet problem in AdS spacetimes with a view to 
understanding the boundary CFT interpretation. We define the problem as bulk Ein- 
stein's equations with Dirichlet boundary conditions on fixed timelike cut-off hyper- 
surface. Using the fluid/gravity correspondence, we argue that one can determine 
non-linear solutions to this problem in the long wavelength regime. On the boundary 
we find a conformal fiuid with Dirichlet constitutive relations, viz., the fiuid propa- 
gates on a 'dynamical' background metric which depends on the local fluid velocities 
and temperature. This boundary fluid can be re-expressed as an emergent hypersur- 
face fluid which is non-conformal but has the same value of the shear viscosity as the 
boundary fluid. The hypcrsurface dynamics arises as a collective effect, wherein effects 
of the background are transmuted into the fluid degrees of freedom. Furthermore, we 
demonstrate that this collective fluid is forced to be non-relativistic below a critical 
cut-off radius in AdS to avoid acausal sound propagation with respect to the hyper- 
surface metric. We further go on to show how one can use this set-up to embed the 
recent constructions of flat spacetime duals to non-relativistic fluid dynamics into the 
AdS/CFT correspondence, arguing that a version of the membrane paradigm arises 
naturally when the boundary fluid lives on a background Galilean manifold. 
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1 Introduction 

The AdS/CFT correspondence [1] which postulates a remarkable duality between large N 
quantum field theories and gravitational dynamics, provides a useful theoretical laboratory 
to address questions underlying the dynamics of these systems. Not only has it proven useful 
to obtain quantitative information about the dynamics of strongly coupled field theories, but 
it also provides a unique perspective into the geometrization of field theoretic concepts. 

Since the early days of the AdS/CFT correspondence it has been known that the radial 
direction of the bulk spacetime encodes in some sense the energy scale of the dual field theory 
[2]. While the nature of this map is not terribly precise outside of the simple example of 
pure AdS geometry (dual to the vacuum state of the field theory), it nevertheless provides 
valuable intuition about certain basic aspects of effective field theory dynamics [3, 4], and 
has led to the idea of the holographic renormalisation group [5], which relates the radial 
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'evolution' in AdS to RG flows in field theories. More recently this idea has been exploited 
to geometrize Wilson's concept of integrating out momentum shells to generate field theory 
effective actions, in terms of integrating out regions of the bulk geometry which in turn lead 
to effective multi-trace boundary conditions on the cut-off surface, a fixed radial slice (in 
some preferred foliation) in AdS [6, 7]. One of the key features of this holographic Wilsonian 
approach was the emergence of multi-trace deformations of the field theory even in the planar 
limit, consistent with field theory expectation. 

A natural question in this context is what does this RG flow mean for the gravitational 
equations of motion? More precisely, consider the problem of integrating out radial geo- 
metric shells in Einstein gravity with negative cosmological constant (which is a consistent 
truncation of string theory /super gravity). One anticipates based on the standard dictio- 
nary which relates the bulk metric to the boundary energy momentum tensor to obtain a 
scale dependent effective action for the energy momentum tensor, containing arbitrarily high 
multi-traces of the stress tensor. The reason for the generation of these multi-traces is clear, 
once one factors in the intrinsic non-linearity of gravity. The basic equations in this context 
are of course easy to write down; as explained in [6, 7] the flow is driven by the radial ADM 
Hamiltonian and one can in principle solve the resulting Hamilton- Jacobi like equation for 
the effective action on the cut-off hypersurface. Despite the conceptual simplicity of the 
formulation of Wilsonian RG in terms of geometric effective actions, the point still remains 
that gravity's intrinsic non-linearity makes explicit solutions hard to come by. 

One can ask whether there is a tractable sector of the gravitational flow equations which 
leads to new insight. A natural avenue for exploration is suggested by the long- wavelength 
regime where we restrict attention to fluctuations of low frequency in the field theory direc- 
tions. As evidenced by the fluid/gravity correspondence [8] there is an essential simplification 
in this regime; bulk Einstein's equations can be explicitly solved order by order in a long- 
wavelength expansion along the boundary.^ As such one should be able to use this framework 
in conjunction with the fluid dynamical expansion to derive an effective action for the low 
frequency degrees of freedom which live on a cut-off surface in the interior of the AdS space- 
time.^ Rather than tackle this problem directly we will take a slightly different tack in this 
paper, one which we believe clarifies some aspects of evolution in the radial direction and its 
possible connection to RG flows. One of our main conclusions will be that imposing rigid 
cut-offs in AdS is more naturally viewed in terms of perturbing the GET by some non-local 
deformation or equivalently by introducing explicit state-dependent sources in the boundary 
theory. 

A second motivation is the recent work [15, 16] which derives an explicit map between 

^See [9, 10] for reviews which describe developments in this area and [11, 12, 13] for generahzations that 
will be of great use in the following. 

^It is natural to expect that such a construction might provide a holographic derivation of the hydrody- 
namic deconstruction described in [14]. 
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solutions of vacuum Einstein equations (with no cosmological constant) and those of incom- 
pressible Navier-Stokes equations, thereby making direct contact with some of the ideas of 
the black hole membrane paradigm in asymptotically flat spacetime [17, 18]. This problem, 
which has been further generalized in [19], is the zero cosmological constant analog of the 
problem we consider (see also [20, 21] for another approach and [22] for related work). The 
idea is to consider a fixed timelike hypersurface with Dirichlet data enforcing a fiat metric on 
the slice. Given these boundary conditions one wants to solve vacuum Einstein equations so 
as to obtain a solution which has a regular future horizon.'^ By explicit construction which 
involves long wavelength fiuctuations around fiat space in a Rindler patch the authors of 
[16, 19] construct solutions to vacuum Einstein's equations order by order in a perturba- 
tion expansion in gradients along the hypersurface directions. The resulting geometry has a 
regular Rindler horizon, and one obtains a regular solution to Einstein's equations contin- 
gent on the fact that dynamics of the induced stress tensor on the hypersurface satisfies the 
incompressible Navier-Stokes equations. 

While this development is fascinating, one is hampered from a first principles under- 
standing of the physics from a holographic viewpoint, owing to the rather poorly understood 
concepts of fiat space holography. Moreover, given the connection between fiuid dynam- 
ics (albeit relativistic and conformal) and Einstein's equations with negative cosmological 
constant as described by the aforementioned fiuid/gravity correspondence [8] (and its non- 
relativistic extension in [13] ), it is interesting to ask whether the construction in [16] can 
be obtained as a limit of the fiuid/gravity map. If this is possible, one can then look for the 
field theoretic interpretation of the fiat space problem. 

Motivated by these issues, we consider a region of the AdS spacetime bounded by a 
timelike hypersurface E^i at some radial position, say r = r^i in the supergravity limit of 
AdS/CFT. We are interested in solving for the bulk dynamics where will give ourselves the 
freedom to specify boundary conditions on S/j. The second boundary condition (which is 
necessary to zero-in onto a unique solution) will be specified by demanding regularity in 
the interior of the spacetime. We have schematically depicted the set-up in Fig. 1. In the 
large N limit, the specification of the problem thus is tantamount to solving classical partial 
differential equations (PDEs) in AdS with a Dirichlet boundary condition imposed on various 
fields at the hypersurface S^). The question we would like to know the answer to is simply: 
'What is the problem that we are solving in the CFT language?" 

As we have reviewed above, various results exist in literature that suggest that solving 
such a Dirichlet problem is analogous to some kind of RG from the CFT point of view. 
Despite the strongly suggestive nature of this holographic RG point of view, it is also not 
very clear what kind of an RG is one speaking of within a CFT. A-priori, for one, it does 

■^In the asymptotically flat spacetime one also has to specify initial data for radial evolution on the past 
null infinity . The boundary conditions chosen in [16] are such that no disturbance propagates into the 
bulk spacetime from J^^ . 
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Fig. 1: Schematic representation of the Dirichlet problem we consider in this paper. The Dirichlet 
surface is taken to be at some value r = ru where we impose boundary conditions on the 
fields. The solutions will further be constrained by requiring that they be regular on any 
putative horizon ^K"*" (shown in the figure) or the origin. The question we are after is what is 
the boundary image of this Dirichlet data? 

not seem like the RG flow tliat arises from cutting off a CFT a la Wilson is the correct way 
to dualize the Dirichlet problem. Hence, our question - what is the CFT dual of a bulk 
Dirichlet problem?^ 

As a warm-up we first consider the bulk Dirichlet problem for linear PDEs, using the 
simple setting of a Klein-Gordon field propagating in a cut-off AdS spacetime. In this case it 
is not hard to see that one is deforming the field theory by a non-local double-trace operator, 
whose precise form, we argue, can be extracted by suitable convolution of appropriate bulk- 
to-boundary propagators. 

We then turn to the issue of setting up the problem in a gravitational setting, outlining 
it in general before moving on to the tractable setting of the fluid/gravity regime. In the 
long wavelength regime we will argue that the bulk Dirichlet problem reduces to a particular 
forcing of the fluid on the boundary of the asymptotically AdS spacetime. The fluid/gravity 
correspondence was generalized to fluids propagating on curved backgrounds with slowly 
varying curvatures in [11] and the most general solutions which will prove to be of interest 

^To avoid confusion we wish to emphasize that we will refer to the bulk Dirichlet problem as one defined 
in the preceding paragraph. While this reduces to the standard Dirichlet problem on the boundary of AdS 
when we take the surface S/j to infinity, we will soon see that the bulk Dirichlet problem induces different 
boundary conditions at infinity when the surface Tijj is retained at a finite position and it therefore pays to 
maintain the distinction. 
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to us were presented in [12]. Using these results it transpires that we can immediately write 
down the solution to the bulk Dirichlet problem in the long wavelength regime. 

The logic is the following: we wish to prescribe on the hypersurface r = a. Lorentzian 
metric which we denote as g^^. This is arbitrary subject to the requirement that its curva- 
tures be slowly varying so that we can treat it with the fluid/gravity perturbation scheme. 
We then solve Einstein's equations demanding regularity in the interior of the spacetime. 
Using standard intuition from the AdS / CFT correspondence it can be argued that the seed 
geometry which we need to set up the gradient expansion should simply be a black hole ge- 
ometry which has a regular future event horizon, which furthermore satisfies the prescribed 
Dirichlet boundary condition.^ It is not hard to see that such a seed solution is obtained by 
simply performing a coordinate transformation of the well known planar Schwarzschild-AdS 
black hole. 

But this is precisely the set-up of [11, 12], the only difference being the fact that in these 
works the Dirichlet data is imposed on the boundary at r = oo. Let's call this boundary 
metric Qf^^, which is also by definition slowly varying etc.. The solution to the asymp- 
totic Dirichlet problem is characterized by the boundary metric Qf^u, a distinguished velocity 
field u'^ (which is unit normalized) and a scalar function b (determining the temperature or 
equivalently the local energy density). Let us denote these variables collectively as X. The 
boundary Brown- York stress tensor (up to counter-terms) takes the fluid dynamical form 
and is built out of the data contained in X. 

Now given the space of solutions to the asymptotic Dirichlet problem, we can reparametrize 
that space of solutions appropriately to obtain the solutions of the new bulk Dirichlet prob- 
lem. The only condition we have to satisfy is that the induced metric on S^) in the solutions 
obtained this way^ be equal to g^j^i,. Furthermore, we can extract the stress tensor on Ej)^. 
We will argue that there is a corresponding velocity field (normalized with respect to the 
hypersurface metric) and a scalar function (the hypersurface temperature), which parame- 
terize the stress tensor of the hypersurface, which not surprisingly takes the fluid dynamical 
form. The main novelty is that the stress tensor does not however correspond to that of a 
conformal fluid. The introduction of an explicit scale by way of the Dirichlet surface's loca- 
tion engenders a non- vanishing trace, which curiously evolves in a highly suggestive manner 

^The implicit idea behind the uniqueness of the seed geometry here is based on the notion that equihbrium 
dynamics in the field theory at finite temperature is governed by a black hole, i.e., we are always in the 
'deconfined phase' in the field theory on the boundary. 

^Note that in order to get something interesting, it is important that we have a solution space allowing 
for non-trivial metrics at infinity (which we do thanks to [12]). We are not simply slicing a single spacetime 
but working with a family of geometries characterized by the boundary data X which is being adjusted so 
as to agree with the Dirichlet boundary conditions. 

^There is an issue of counter-terms that one can use when working at finite radial coordinate in an 
asymptotically AdS spacetime. We will take the conservative view that the relevant counter-terms are the 
same as those necessary asymptotically. 
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under change of cut-off surface position, see (3.26). 

Calling the totality of the data on the hypersurface X we further show that within the 
gradient expansion there is a one-to-one correspondence between the hypersurface data and 
the boundary data; X X is bijective. This then has the advantage that we can 

immediately understand the boundary dual of the bulk Dirichlet problem as a conformal 
fluid which is placed on a^ 'dynamical background' whose metric depends on the same set of 
variables that characterizes the fluid itself (in addition to the prescribed hypersurface metric). 
So from the boundary viewpoint there is a complete mixing between intrinsic and extrinsic 
data, which is the long-wavelength non-linear analog of the double trace deformation seen 
for the scalar toy model. Moreover, this solution allows us to see that the dynamics of 
the fluid on the Dirichlet surface, as given by the conservation equation on S^, 'emerges' 
as collective dynamics of the boundary CFT. In particular, the boundary fluid lives on a 
'dynamical background', and the effects of the background can be suitably subsumed into the 
fluid description. This suggests that the correct way to think about the hypersurface physics 
is in terms of a 'dressed fluid' living on an inert geometry. Thus, the effective description of 
a fluid on this dynamical background is geometrically encapsulated in terms of the Dirichlet 
hypersurface dynamics. 

Examining the resulting dynamics on Ed we find that the hypersurface or effective fluid 
suffers from a possible pathology for smaller than some critical value r£,^snd- At r£,^snd 
the sound mode of the effective fluid starts to propagate outside the inert background E^'s 
light-cone. We suggest that in the CFT, this effect is due to the extreme forcing of the 
fluid on the boundary by the 'dynamical' metric, and moreover propose that one can obtain 
sensible dynamics by projecting out the sound mode. This involves looking at the fluid at a 
scaling limit and this can be formalized as taking the incompressible non-relativistic limit of 
the fluid on the hypersurface in a manner entirely analogous to the scaling limit described 
in for generic relativistic fluids in [13, 23].^ 

Having understood the Dirichlet problem for generic E/) away from the horizon 
we then proceed to push this surface deeper into the spacetime and ask what happens 
as we approach the horizon. In this regime E/) dynamics continues to be described by 
incompressible Navier-Stokes equations in the limit, though with some slight differences 
from the BMW limit mentioned above. Zooming in onto the region between E/) and the 

®We use the word dynamical to characterize the boundary metric in the following sense: the metric on 
the boundary depends on the dynamical degrees of freedom of the system, viz., and T, as in (4.1). This 
is a pre-specified constitutive relation for the boundary metric in terms of the fluid variables, not unlike the 
constitutive relation for the energy momentum tensor. We will call such a constitutive relation coming from 
bulk Dirichlet problem as the Dirichlet constitutive relation. The fluid at the boundary, therefore, sees a 
dynamic metric background but no new degrees of freedom are introduced and hence for example, there are 
no boundary Einstein's equations that need to be solved. 

^We shall refer to scaling limit as the BMW limit after the authors of [13], despite the fact that we are 
focussing on sub-sonic excitations and that it has been well documented in classic textbooks, cf., [24]. 
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horizon, we provide an embedding of the construction of [16, 19] into the fluid/gravity 
correspondence [8]. Further, we demonstrate that in this hmit both the bulk metric in the 
region between and the boundary, and the boundary metric degenerate from metrics on 
a Lorentzian manifold to Newton-Cartan like structures. This raises interesting questions 
about the natural emergence of the Galilean structures in the AdS/CFT correspondence 
which we postpone for future work. 

The plan of this paper is as follows: In §2 we first address the Dirichlet problem for a 
scalar field propagating in AdS^+i using this linear problem to build intuition. In §3 we pose 
the bulk Dirichlet problem for gravity in AdS spacetime and solve it in the long wavelength 
approximation borrowing heavily on the results from the fluid/gravity correspondence. The 
remainder of the paper is then devoted to understanding the physics of our construction 
in various regimes: §4 demonstrates how the Dirichlet surface dynamics, as governed by 
the conservation equation, arises from the boundary physics. Aided by this we argue that 
the Dirichlet dynamics is probably pathological past a critical radius and propose a non- 
relativistic scaling of the resulting fluid a la BMW in §5 to cure this possible pathology. 
Finally, in §6 we study the near-horizon Dirichlet problem and make contact with the recent 
work on the flat space Dirichlet problem (and its connection with Navier-Stokes equations). 
We end with a discussion in §7. Various appendices contain useful technical results. In 
particular, to aid the reader we provide a comprehensive glossary of our conventions and 
key formulae in Appendix A. This is followed by a complete 'Dirichlet dictionary' relating 
hypersurface variables to boundary variables in Appendix B for ready reference. 

Note added: While this work being completed we received [25] which has partial overlap with 
the results presented in §3. These authors also attempt to solve for bulk geometries with 
prescribed boundary conditions on S/j in the long wavelength regime and interpret their 
results in terms of a RG flow of fluid dynamics. 

Note added in v2: In the flrst version of the paper the non-relativistic metrics quoted in §5 
and in Appendix D were incorrect; the metrics as presented do not solve the bulk Einstein's 
equations to the desired order. These are now corrected in the current version. However, 
the full set of terms that we need to include in order to see the Naiver-Stokes equation on 
the boundary is quite large. Hence in the main text we only report the results for the case 
where the non-relativistic fluid moves on a Ricci flat spatial manifold in §5 and present the 
general results in a new appendix Appendix E. We note that the results of §D also correct 
the expressions originally derived in [13]. 

2 Dirichlet problem for probe fields 

To set the stage for the discussion let us consider setting up the bulk Dirichlet problem 
for linear PDEs in an asymptotically AdS^+i spacetime. As a canonical example we will 
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consider the dynamics of a probe scalar field $(r, x'^) of mass m. Generalizations to other 
linear wave equations such as the free Maxwell equation are straightforward. 

We will let this scalar field propagate on a background asymptotically AdS^+i geometry 
with spatio-temporal translational symmetries so that the background metric can be brought 
to the form 

ds^ = g^,{r) dx^" dx" + dr^ (2.1) 

The boundary of the spacetime is at r — )■ oo and we will assume that the boundary metric is 
the Minkowski metric on M*^"^'^ for simplicity, so that asymptotically g^j^^, — )■ rj^i, and Qrr 1 
(as r — 7- oo). 

The dynamical equation of motion for the scalar is the free Klein-Gordon equation which 
can be written as an ODE in the radial direction for the Fourier modes ^k{f) of $(r, x^) 

^r,xn = J ^''■^'^^W, (2.2) 

and takes the form 

, ^ dr (v/=^9,<l>fc(r))-((7^,Pr + m2)<|., = (2.3) 

v-ggrr 

As a second order equation we need to specify two boundary conditions. We are going 
to restrict attention to the finite part of the geometry as illustrated in Fig. 1 and impose 
Dirichlet boundary conditions for the field at some hypersurface E^^ at r = r/j. The second 
boundary condition in general can take the form of a regularity boundary condition in the 
interior of the spacetime. If we were working in a spacetime with a horizon this would 
demand that the mode functions of interest are purely ingoing at the future horizon. 

The question we wish to pose is the following: usually in an asymptotically AdS(i+i 
spacetime we know that the solution to the scalar wave equation above has two linearly 
independent modes with power-law fall-off characterized by the source and vev of the 
dual boundary operator 0$ (which we recall is a conformal primary). We wish to ask what 
is the characterization of the boundary data as a functional of the Dirichlet hypersurface 
data. In this simple linear problem it is easy to see that there is a one-one map between the 
two sets of data. Essentially we are asking how to tune and so that the value of the 
scalar field on the Dirichlet hypersurface a.t r = takes on its given value. 

While it is possible to derive a formal answer to the above question, it is useful to first 
visit the simple setting of pure AdS^+i spacetime where we have the luxury of being able to 
solve the scalar wave-equation explicitly to see an explicit answer to the question. 

2.1 Probe scalar in AdS^+i 

We specialize our consideration to the pure AdS^+i geometry where (7^^ = rj^i, and grr = 1 and 
one has enhanced Lorentz symmetry on the constant r slices. The wave equation simplifies 
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T^Tr ('■^« |*.-W) + . "=)1>. = (2.4) 

This is well known to have solutions in terms of Bessel functions, but we will proceed to 
examine the behavior in a gradient expansion to set the stage for the real problem of interest 
later. 



2.1.1 The k = case 

The translationally invariant solution of the massive Klein-Gordon equation (2.4) in the bulk 



is^o 

*W = ^^ + r^-% (2.5) 

where the dual primary has a scaling dimension A obeying A(A — ci) = along with a 
source and a normalized vev^^ defined via 

J,^K-^$(r)]_ (2.6) 

^ ^ (^d-A^)]^^^ ^ ^_^A ^^Q^^ _ _ (2.7) 

with 

z.^A-- = ^^+m2 (2.8) 

for convenience. For simplicity, we will assume A > | and choose m such that z/ ^ Z to 
avoid complications with logarithms. Extension to A G [| — 1, |] with the lower end of 
the interval saturating the unitary bound is possible with the added complication of taking 
proper account of the necessary boundary terms. 

We will begin by rewriting this solution in terms of the quantities on the Dirichlet surface 
which we denote with a hat to distinguish them from the boundary data: 

J,^[r^-^$(r)]_,^ = J, + ^^, (2.9) 

^ ^ (^d-A^)]^^^^ ^ ^ _ ^2.10) 

Since the transformation between the data on the boundary {J^^ip} and that on the 
hypersurface {J^, 0} is linear it is a simple matter to write the bulk solution in terms of the 
hypersurface variables. One simply has 



$(^) = + ^A-d J _ H_ .2.11) 

^ ^ (2z/) [ ^ 



^''We will for the moment refrain from imposing any IR boundary condition so as to be able to see the 
general structure. After all in pure AdS imposing regularity at the Poincare horizon would kill the vev which 
has to vanish in the vacuum. 

^^Note that the vev is (IGtt G^+i)^^ (t> where Gd+i is the gravitational constant in AdS^+i. 
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This is the answer we seek and all that remains is to interpret this result. 

It is now easy to notice that the imposition of the Dirichlet condition on a hypersurface 
inside the bulk is equivalent to making the boundary source a specific function of the vev. 
From (2.11) we can read off the specific deformation of the boundary CFT action to be given 
by 

which happens to be an irrelevant double-trace deformation [26, 27] of the boundary CFT. 
Hence, at least in this simple setup the dual of the Dirichlet problem is to make the source 
of a primary 0$ a particular joint function of the vev of the primary in the given state and 
another fixed (state-independent) auxiliary source. 



2.1.2 The k ^ case : Derivative expansion up to k"^ 

Having seen the result for the translationally invariant case = 0, we now proceed with 
k ^ 0. It is well known that general solution to the wave equation (2.4) is given in terms of 
Bessel functions which we parameterize as^^ 

= ^ X ^^.''W^) + -"-'M. X ^I^KAk/r) (2.13) 

Note that our previous result for A; = follows from just keeping the leading terms in the 
expansions 

r(^) ^ ^ m j\ ^ T{u) j\ 

' ' (2.14) 

2x- ^ ^(2z/ + 2j)r(z. + j) j! 

For a general k, we can repeat the analysis of the previous section. While this can be 
done generally at all orders in k with some work, for simplicity we will resort to derivative 
expansion keeping terms upto order fc^. Not only will this allow us to see some of the 
structures emerging explicitly, but it also sets the stage for our gravitational computation in 
later sections. 



^^We are working here with space-Hke momenta having > 0. Results for time-hke momenta follow from 
replacing fc by i fc and using the Bessel relations 

2{ixY IT 2x^ 

K,{2ix) = ^^Y,{2x) + i sin [{-i^)^] T{-v + l)x'' M^x) 



2{ixY 2x'' 



11 



Using the expansion (2.14) above, we have 

2u P 



(2z/)r^ 



1 + 



(2z/ + 2)2 2r2 

2v 



+ . . . I +r 

1 

1 



f2z/ + 2)2r2 



{J(f>)k 

+ 



{2u - 2) 2r^ 



+ 



(2.15) 



with the eUipses representing order fc^ terms and higher. 

The source at the intermediate surface is as before easily determined 

= {J^)k I 1 



1 



+ 



{2u-2)2rl 
Jk I ^ 2z/ 



fc2 T(-i^) 
+ . . . + 



{2u) rl^ 



2r 



D 



2v 



{2u + 2)2rl 



(2z/ + 2)2 2r2, 



(2.16) 



while the normalized vev of the primary to this order in derivative expansion can be deter- 
mined after subtracting an appropriate counter-term as^^ 



2z/-2^ 



1 + 



2z/-2 P 



(2z/)2 2r 



D 



I 4(J^)fe 



(2i/-2) r(z/ + 2) V2 



2u 



(2.17) 



To solve the Dirichlet problem,we need to solve for 0^, (J^)^ in terms of the hatted 
variables from (2.16) and (2.17) which can be inverted to get 



(Pk_ 
2) 



(2z/)2 2rl 

1 e 

{2u-2)2^ 



+ 



D (2z/) rj^ 
T(-u) ( k 



2v 



„2 
D 



D(2z/-2) r(z/ + 2) V2 



+ . . . 



2v 



(2z/ + 2)2 2r2 

1 A;2 
{2v + 2)2rl 



+ 



(2.18) 



^■^If we keep fc* terms and higher, one needs to subtract appropriate counter-terms at that order. These 
counter-terms are determined by requiring that 4>k is finite as r^i — ^ oo. The expUcit expressions for counter- 
terms to any required order can be determined using the expansions in (2.14) (see for e.g., [28]). 
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where the momentum dependent coefficient T) is 



(2z/)2(2z/-2)2r|5 r{u) \2rD) \ u^{u - 1) 
(2z/ + 1)^ - 4(2;/ + 2)2 + 7 ^ 



(2z/)2(2z/ + 2)3 r2 



(2.19) 



As we saw in the k = case, we have yet again determined a state dependent source on 
the boundary for the primary operator 0$. The key feature to note from the above analysis, 
is that the expression for the boundary source is non-analytic in k and hence non-local 
when Fourier-transformed back to position space. Hence, we see that in general we have 
a map between the non-local double trace deformation on the boundary and the Dirichlet 
data on S/) (similar non-local double-trace deformations were explored earlier in [29]). 



2.2 A general proposal for linear systems 

From the analysis of the free scalar wave equation in AdS^+i the picture is rather clear. In 
the CFT, in general one can make the source a non-local functional of the vev of the primary 
operator. Usually such a function can be fed into the holographic dictionary via a 'state- 
dependent' boundary condition, which whilst somewhat unnatural from a field theory is a 
perfectly sensible boundary condition to consider. For some special classes of functionals, 
this state-dependent boundary condition has a very simple bulk interpretation as a Dirichlet 
boundary condition imposed on an intermediate surface, implying that we can trade the 
non-locality of the boundary sources into local behavior at some lower radius. 

We just have one further question to answer before we declare victory: how do we in 
practice determine this special set of sources in various holographic setups? For the general 
backgrounds we can formally write the solution to the wave equations in terms of integrals 
over the Dirichlet data convolved with suitable 'Dirichlet bulk to boundary propagators', 
3Csource a^d X^cv The former propagates the information contained in to the boundary 
source, while the latter allows determination of the contribution from the vev on S/j, i.e., 
formally 

J<f>{x) = j d'^x' ^^XsourccirD,x;x') J^{x') + %^cv{rD,x;x') (j){x')^ (2.20) 

Note that implicit in our definition of these Dirichlet bulk to boundary propagators is the 
information of the boundary condition in the interior of the geometry and the necessary 
counter-terms. While it is possible to work this out in more specific geometries, such as a 
Schwarzschild-AdSrf+i spacetime to see the interplay of these IR boundary conditions, we 
will leave this toy problem for now, and proceed to analyze the more interesting case of 
gravitational dynamics in AdS^+i wherein we do have to face-up with non-linearities of the 
equations of motion. 

^^In fact, even for pure AdS^+i it is interesting to ask what the deformation on the boundary is when we 
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Before proceeding to the gravitational setting, however, let us make a few pertinent 
observations relevant to the motivation mentioned at the beginning of §1. The result we have 
obtained is quite intuitive; demanding that our fields take on the desired value at S^) entails 
a linear relation between the two pieces of data at infinity, thereby leading to the observation 
about the source depending on the vev. We also see that despite some superficial resemblance 
to the Wilsonian RG flow where too one encounters multi-trace operators there is a crucial 
distinction in the physics. In the formulation of [6, 7] one finds that for fixed asymptotic 
data, upon integrating out the region of the geometry between the boundary and a cut-off 
surface (which we can for simplicity take to be S^i for the sake of discussion) one obtains an 
effective action for a cut-off field theory living on S/j with scale dependent sources. These 
are irrelevant double traces (which are the only terms generated in a Gaussian theory which 
the linear models under discussion are), and one obtains the /3-functions for the double trace 
couplings along the flow. In the present context however what we have is a situation wherein 
we are forced to engineer a specific double trace deformation on the boundary so as to ensure 
that we satisfy the Dirichlet boundary conditions on T^d. This is conceptually different from 
usual notions of RG, where one does not conventionally consider state dependent boundary 
conditions in the UV. However, there is a sense in which renormalisation of sources takes 
place which will become quite clear when we look at the gravitational problem. 

3 The Dirichlet problem for gravity 

Having understood the boundary meaning of the Dirichlet problem for probe fields in a fixed 
background, we now turn to the situation where we consider dynamical gravity in the bulk. 
While we could consider other matter degrees of freedom in the bulk whose backreaction 
we now have to take into account, we choose for simplicity to restrict attention to the 
dynamics in the pure gravity sector which, as is well known, is a consistent truncation of 
the supergravity equations of motion. From a field theory perspective, we are going to work 
in the planar limit and focus on the dynamics of a single operator, the stress tensor and its 
source, the CFT metric g^j^^,. 

3.1 Setting up the general Dirichlet problem 

First of all one should ask what does it mean to consider the Dirichlet problem at a fixed 
hypersurface in the bulk when gravity is dynamical. We will take the view that the location 
of the hypersurface E/) is specified by a scalar function on the bulk manifold M^+i. We want 
to determine the metric on this spacetime by solving the dynamical equations of motion 
subject to the boundary condition on the prescribed hypersurface. To wit, we demand that 

move S^) close to the Poincare horizon; in this Umit it seems natural to expect that the asymptotically we 
obtain a a Neumann boundary condition. We thank Don Marolf for emphasizing this to us. 
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be endowed with a Lorentzian metric 9mn which solves Einstein's equations with a 
negative cosmological constant. The equations of motion are (in units where RAds = 1)^^ 

Emn = ^MN + d 9mn = (3.1) 

We will adapt coordinates = {r, x^} to the hypersurface T,d and take this distinguished 
surface to be at r = r^j with intrinsic coordinates x'^.^^ We impose the boundary condition 

Smn dX^' rf^^^l,^,^ = rl g^,{x) dx^ dx" (3.2) 

where g^Mv^x) is the Dirichlet data we wish to specify. To complete the specification of 
the problem we should impose some boundary condition in the interior of the spacetime, 
which we will canonically take to be a regularity condition. The scaling by r|, above whilst 
unconventional from the bulk perspective, is more natural in the AdS/CFT context for it 
makes it easy to compare with the case where we push the hypersurface to the boundary. 

The general problem as stated above is quite hard. For one it is not clear that for generic 
choices of Dirichlet data one obtains a solution compatible with regularity in the interior 
of the spacetime. While a local solution in an open neighbourhood of S/j can presumably 
be obtained by adapting a Fefferman- Graham like expansion, one is unlikely to be able to 
gain much insight using this procedure. Moreover, given that the map between Ejj and 
the boundary is expected to be non-local (borrowing intuition from the linear problem) one 
wonders whether there are causality issues as well. In particular, for generic state dependent 
boundary conditions causality is murky - does the source adjust itself acausally to obtain 
the appropriate response? Likewise on there is a concern that signals can propagate 
outside the light-cone of the metric g^^ (which they could for instance do through the bulk); 
does this imply a corresponding pathology for the boundary physics as well?^^ In short, the 
well-posedness of the Dirichlet problem is a-priori unclear. Nevertheless, we will ignore all 
these subtleties for now and forge ahead. 

After solving Einstein's equations with the boundary conditions we have set-up, we can 
extract the Brown- York stress tensor on S/j, denoted T^^, using the standard set of boundary 
counter-terms [30, 31] 

= -^^^ {k^.u -Kg,, + {d- 1) g,. + ...). (3.3) 

rjy K^y is the extrinsic curvature of the hypersurface S £, and K is its trace, defined as usual 
in terms of the normal to the surface. 

^^We will use upper-case Latin indices for the bulk spacetime indices, reserving lower-case Greek indices 
for hypersurface or boundary indices. See Appendix A for a list of conventions. 

^^It might be useful to view this scalar function as physically being specified either as the level set of the 
red-shift factor, or by introducing a dynamical scalar field. 

^^We thank Veronika Hubeny for extensive discussions on these issues. 
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Note that we have written the answer in terms of the intrinsic metric on the hypersurface 
which is related to the induced metric from the bulk up to a rescaling by r|). The homoge- 
neous scaling of the stress tensor under allows us to fix an overall rr, dependent pre-factor. 
The hypersurface stress tensor (3.3) is of course covariantly conserved (the Gauss-Codacci 
constraint), i.e., 

V^T^ = 0. (3.4) 

Given {g^^,T^i^} we can ask what are the corresponding boundary conditions on the 
boundary that lead to the same geometry. Based on our scalar problem we can conclude 
that the boundary source g^^ and stress tensor T^^, are in general non-local functional of 
the Dirichlet data. We would like to characterize the map between these two sets of data^® 

<Pd : {g^lu,f^u} -> {gf,u,T^u} ■ (3.5) 

While this problem is in general difficult, there is one context in which we can not only 
solve for the boundary data in terms of the hypersurface variables, but we can also investi- 
gate the issues raised above in precise terms. This is the long-wavelength hydrodynamical 
regime along the hypersurface as in the fiuid/gravity correspondence [8, 11, 12], wherein 
gravitational duals to arbitrary fiuid fiows on the boundary were constructed order by order 
in a gradient expansion. The reason this is possible, as explained in these works, is that 
the bulk spacetime in this long-wavelength regime is well approximated 'tubewise' by the 
boosted planar Schwarzschild-AdS^+i solution, see Fig. 2. As a result one finds that Ein- 
stein's become ultra-local in the directions leading one to effective ODEs to determine 
the radial profiles. The tubes in question are centered around radially ingoing null geodesies, 
which can be used to translate information from any bulk hypersurface of interest to the 
boundary (see [32, 12] for a discussion of the causal structure). Given this, it is actually 
easy to solve the problem of finding the map (fo and we now describe the construction in 
the rest of this section. 

3.2 Dirichlet problem and the Fluid/Gravity correspondence 

To set the stage for our discussion, let us recall that the fiuid/gravity map constructs in a 
gradient expansion, regular solutions to the bulk Einstein's equations (3.1) which are dual to 
arbitrary fiuid fiows on the boundary of the asymptotically AdS^+i spacetime. For a bound- 
ary metric g^v{x) which is slowing varying, the boundary stress tensor in this context is a 

At this point it is worthwhile to get a technical point out of the way. We will have two metric structures 
in the story henceforth, the hypersurface metric 5^1/ and a boundary metric 5^1^. To avoid confusion, we will 
write all equations intrinsic to the hypersurface or to the boundary consistent with the respective metric 
structures. In practice this simply means that one raises/lowers indices of the equations with respect to the 
appropriate metric. These have to be handled with care, but by judicious use of the two metrics one can 
relate other components if required. 
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Fig. 2: Schematic representation of the gravitational Dirichlet problem in the fluid/gravity regime. 

The causal structure of the fluid/gravity spacetimes is illustrated emphasizing the tubewise 
approximation; in each tube the geometry resembles that of a uniformly boosted Schwarzschild- 
AdS(i+i black hole. Suitable choices of the Dirichlet surface allow us to find the map between 
the boundary data X and the Dirichlet hypersurface data X within each tube, rendering the 
problem tractable. 

not an arbitrary symmetric traceless two tensor, but constrained to take the hydrodynamical 
form. It is parameterized by d independent parameters - a velocity field (unit normal- 

ized so that g^y = —1) and a scalar function b{x) which parameterizes the temperature. 
The bulk metric Smn is determined in terms of the data X = {g^u{x),u^{x),b{x)}. 

We wish to implement the same procedure, but starting with analogous data X = 
{(jfiuix), Uf_i{x), b{x)} on the Dirichlet hypersurface T^d. But given the ultra-locality inherent 
in the long-wavelength regime and the fact that [12] have solved the problem for arbitrary 
boundary metrics (corresponding to fluids on arbitrary slowly varying curved backgrounds), 
we don't need to solve any equations. The solution space of the bulk Dirichlet problem 
coincides with the solution space found in [12] and the problem at hand is readily solved 
by slicing this solution space appropriately. With this aim, we now review the solutions 
constructed in [12]. 
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3.2.1 Review of fluid/gravity 

The general solutions of the bulk equations of motion (3.1) in the fluid/gravity regime take 
the form 



^MNdX^^ dX^ = -2\x^{x)dx^ (dr + rQJ^(r,a;) dx"") + <5^^{r,x) dx^" dx"" , (3.6) 



where the fields 03^ and S5^jy are functions of r and x^ which admit an expansion in the 
boundary derivatives and are known to second order in the gradients. For our purposes it 
will suffice to consider the first order metric where 



u 



2J;. = ^. + -/(Hw. 







P^, + 2bF{br) a 



flU 



(3.7) 
(3.8) 



with the functions 



1 - 



X 



d ' 



F{x) 



OO yd-l _ 1 

yiV" - 1) 



dy. 



(3.9) 



Afj^ is the Weyl covariant connection introduced in [33] which is expressed in terms of the 
acceleration and the expansion of the velocity 



= u^V\u^ tU„ = a,, — 



Vxu^ 9 

— u,, = a,, ; 

d-1 ^ ^ d-1 



while (j^jy is shear strain rate tensor of m^: 



zzz P " P/3 



d-l 



with P, 



(3.10) 



(3.11) 



So the bulk metric to first order in derivatives explicitly takes the form: 



ds' 



-2 Uf, dx^" {dr + r A^dx") + 



dx^'dx'' + ... (3.12) 



and we have refrained from explicitly denoting the x^ dependence of j£ and the ellipses 
denote second order and higher gradient terms. The corresponding co-metric (the metric on 
the cotangent bundle/the inverse metric) is given by 



f{br) — 



2r9 
d- 1 



dr®dr + 2 [u^ - r-^a^'] 0, dr 



(3.13) 



+ ^-2 [pM- _ 2hF{hr) o^'\ ® 



^^In the fluid/gravity literature one chooses to maintain = to all orders in the gradient expansion 
for simplicity. We will generalize this suitably in the rest of the discussion to simplify our formulae. 
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The stress tensor on the boundary is that of a viscous relativistic fluid: 



T^^y = pgf,u + {e + p)u^Uy- 27] a^,^ + ... (3.14) 
with thermodynamic state variables 

and shear viscosity 

The expression in (3.12) should be thought of as a way to generate solutions of Einstein 
equations when provided with hydrodynamic configurations that solve the ideal fluid equa- 
tions derived form the T^i, above, i.e., when provided with that satisfy V^T'^^ = to 
first order in gradients. Our aim is to reformulate this set of solutions as solutions to the 
bulk Dirichlet problem. 



3.2.2 Dirichlet data from fluid/gravity solutions 

Given the solution (3.12) to the bulk equations of motion, we will begin by simply slicing it 
at a given radial position r = r^, and extract the intrinsic metric on the Dirichlet surface^°. 

The advantage of working with the Weyl covariant form of the bulk metric is that one can 
simultaneously deal with Dirichlet surfaces specified by slowly varying functions r = p{x). 
These can always be brought to the form r = r^) by working in a suitable boundary Weyl- 
frame (local rescaling by a conformal factor log(l — p{x)/rD) will do the trick). The Weyl 
connection A eats p{x) so that 

with this understanding all our formulae hold for arbitrary p{x). 

At fixed r = To the hypersurface metric reads (recalling (3.2)) to first order 

Qfiv = Qiiu + + 26F(&rB) (T^,u - ^M(^yi^) + . . . (3.17) 

While this relation was obtained by slicing a known solution with prescribed asymptotic 
boundary conditions, it has the nice feature of having solved the equations of motion (to the 
desired order) and moreover satisfies the regularity condition in the interior. We will now 
turn the logic around and imagine g^y to be specified at S^) and view the equation above 
as specifying the boundary intrinsic metric in terms of the hypersurface metric. Hence, the 



''We would like to thank Sayantani Bliattacharyya for collaboration on some of the ideas in this section. 
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above equation giving gfj_„ in terms of gf^„ and is the Dirichlet constitutive relation we seek 
- when such a relation is imposed on the boundary data, the intrinsic hypersurface metric 
is automatically fixed to Qf^u- 

To complete the specification of the map (f£, we need to further eliminate the boundary 
velocity field in favour of a vector field on the hypersurface. To do this we need to examine 
the hypersurface stress tensor and parameterize it appropriately. The stress tensor at vd is 
easily obtained from (3.3) to be 

1 2 

T^,u = pQiJiu + ^{e+ p)u^Uy -2ar]a^y^ {e + p)u(^^Au) + . . . (3.18) 

withal 

with T) as given before in (3.16) and we have defined a hypersurface scalar 

d = , ^ = = , ^ (3.20) 

The stress tensor is tantalizingly similar to that of a viscous fiuid, but as yet, we cannot 
interpret 7] as the shear viscosity since we have not expressed T^y, in terms of hypersurface 
variables. 

This is however easy to remedy. Define to be the unit normalized (with respect to g^y 
of course) timelike eigenvector of T^i/. A simple computation shows that 

+ — (3.21) 

a To 

We want to express the hypersurface stress tensor in terms of and its gradients with respect 
to the g^i, compatible connection V^. This can be done by the standard computation of the 
difference of — V^. We outline the calculation in Appendix C and simply quote the 
relevant result here: 

(^riv — CiKTuy , A-i/ — A-u — - ^ , „ — O'u ■ (3.22) 

1 + |(a^ - 1) 

Armed with this data we can write now the stress tensor at r = as 

T^y = pg^^u + {e + p)u^Uy-2ria^y + ... (3.23) 

We see that the result indeed is a stress tensor of a relativistic fluid with energy density e, 
pressure p, given in (3.19) and the same value of shear viscosity as the boundary theory r], 
(3.16). The dynamical content of this system is still the conservation equation 

V^T^. = . (3.24) 

^"'^The expressions for the hypersurface energy density and pressure i and p can be re-written in terms 
of the hypersurface temperature T (see (4.13)) which is the more natural quantity on Ejj. However, it is 
convenient for practical reasons to leave these definitions in terms of b. 
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which follows from realizing that this is the 'momentum constraint' equation for the radial 
slicing of Einstein's equations (or if one prefers the Gauss-Codacci constraint on the hyper- 
surface). Its equation of state is given in (3.19) and we will momentarily evaluate the trace 
of f^, in (3.25). 

Further note that the stress tensor has no contribution associated with the expansion of 
the fluid, i.e., the bulk viscosity vanishes identically on the hypersurface. Nevertheless, the 
fluid is not a conformal fluid, for the trace of the stress tensor is non-vanishing. 

d{d — 1) a 



-e + {d- l)p 



1 a 



(3.25) 



167r Gd+i a + 1 h'^ ' 

This is not entirely surprising for we have introduced an explicit scale ro into the problem, 
and as required for consistency the trace vanishes in the limit r^)— >ooasQ;— t-I. More 
curious is the fact that rate of change of the trace with the S/^'s radial location is simple: 

de 



T% = -ro 



dr 



D 



(3.26) 



The evolution of the trace is highly suggestive for (3.26) can be interpreted as saying the 
the trace is generated by the variation of the local energy density with respect to some 
scale. This kind of a relation probably hints at some kind of non-linear realization of scale 
invariance. Again this is reminiscent of the holographic RG ideas and it would be interesting 
to flesh this out in greater detail. 

Having the notion of the hypersurface velocity fleld we can now proceed to write the 
boundary metric in terms of hypersurface data. Inverting the relation for the velocities to 
obtain (see Appendix C) 



a 



1) 



l + f(d2 



1) 



one can show that 



Qtiy = Qiiu - (a^ - 1) Uy + 



2^2 



1 + i(d2 



1) 



U(,,a 



(3.27) 



2b 

--F(hrD)a^, (3.2^ 
a 



The equations (3.27) and (3.28) together specify the map v^d we seek in the long-wavelength 
regime. Note that the hypersurface and the boundary data are determined by the same 
scalar function h{x), which however enters non-trivially through a in the determination of 
dynamics on S/j. 

Note that the light-cones of g^j^y are enlarged by a factor determined by a relative to that 
of g^u- This is the first signal that there is some interesting interplay between boundary 
causal structures and fixing boundary conditions on E^j; we will address this issue in some 
detail in §4. But first we finish the solution to the Dirichlet problem as stated and write 
down the bulk metric in the long wavelength regime. 
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3.3 Bulk metric in terms of Dirichlet data 



Given the map in (3.27) and (3.28) we are in a position to re-write the bulk metric (3.12) in 
terms of E/j data X alone. Substituting the transformations we can write the final result as 
in (3.6) with 





= 








= i 



a u, 



a 



[l + f(d2-l 



e ^ 

u.. 

d-1 ^ 



f{br) Uf, 



[l + |(d2-l 
P^, + 2hF{hr) a^. 



d 



1 ^ 



P^, + 2bF{br) & 



(3.29) 



and we have defined 



£-l-^^^ = l-f ^/(frr) (3.30) 

2 td j{orr,) 2 ro 



F{hr) = ^{F{hr)-F{hrj,)) = ^ / , ^, dy . (3.31) 



and 

& jbr yiy" 

The factors of are distributed between and 2J,^ by the requirement that the former 
be transverse to u^. 

This bulk metric Smn solves the gravitational Dirichlet problem in the long- wavelength 
regime. We have thus solved the problem posed at the beginning of this section completely in 
terms in this regime aided by the ultra-locality of the gradient expansion. We summarize the 
complete map and the resultant dictionary between CFT variables and the hypersurface 
variables in the Appendix B. As a consistency check note that the results agree when we 
send the surface S/j to the boundary with those derived in [12]; one simply sets vd = oo 
and a = 1. In the remainder of the paper we will explore the relation between the dynamics 
on E/) and that on the boundary, with an aim towards getting better intuition for various 
issues raised in 51 and ^53.1. 



4 Emergence of Dirichlet dynamics in the CFT^ 

Having obtained a solution to the gravitational Dirichlet problem in the long-wavelength 
regime, we now turn to analyze the underlying physics of the system. From the viewpoint of 
the bulk, the dynamics of the system under study has two equivalent descriptions - one in 
terms of hypersurface (hatted) variables V ^T^^ = 0, and another in terms of the un-hatted 
boundary variables V^T^'^ = 0. While the former is more natural in the bulk (since it is 
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Qfj^i, which is fixed in the bulk for our Dirichlet boundary conditions on S^), only the latter 
has a straightforward interpretation in the CFT^. Hence, it is interesting to ask how the 
hypersurface description should be interpreted within the CFT. 

From the CFT point of view, the fiuid is living in a metric background with a Dirichlet 
constitutive relation (rewriting (3.28) in terms of m^) 

/I 29 \ 2 , „ X 

g^lu = Qiiu - [I - ^ + 73 — 7T u^Uu H Ul^^a^) - 2 bF{brD) a^^ + ■ ■ ■ 



gt^u _ ^M^^ _ / 1 _ ^2 _ ^ ^ ^ u^'u" + 2 bF(brD) af"' - — u'^^a"^ + ... 

[d-l)roJ Td 



(4.1) 



where is defined in (3.20).^^ The function a increases with increase in the local temper- 
ature (decreases with increase in the local b). These expressions basically tell us how the 
ambient spacetime background the CFT^ lives on responds to the motion of the CFT fiuid. 

Let us first understand the physical content of the piece in the constitutive relation with 
zero-derivatives. This can mostly be done heuristically which we will do first and then 
confirm it with an explicit calculation. The zero-derivative piece of the boundary metric 
Qfj,^ is made of an inert piece g^^, which does not respond to the fiuid motion along with 
an additional piece proportional to u^Uu- The presence of a term proportional to Uf^Uiy 
means the boundary metric effectively has a correction in its dt"^ piece in the local fiuid 
rest-frame (which is responsible for opening up of the light-cone in the boundary). This kind 
of correction as is well known in general relativity just represents a gravitational potential 
well. Note that this potential well travels along with the fiuid and hence it is tempting to 
think that there is a way to describe the collective packet of fiuid and the local graviton 
cloud that it carries along in terms of a 'dressed' fiuid. 

To guide our intuition, let us draw analogies with another familiar physical situation 
where the background responds to the system locally via this kind of a relation. One analo- 
gous situation is that of a charge carrier moving in a polarizable medium. The polarizability 
of the medium defines the constitutive relation of the medium in exact analogy with the 
constitutive relations for the metric above. We know that in the case of the charge carrier 
moving in a polarizable medium often the polarizability can be taken in to account by shift- 
ing the dispersion of the charge carrier and pretending that this 'dressed' charge carrier is 
essentially moving through an inert medium. 

What we want to argue out in this section is the fact that a similar 'dressing' phenomenon 
happens in the case of the CFT^ fiuid - we want to rewrite the problem of a fiuid with T^'^ 
moving in the 'polarizable' g^i, into the problem where a dressed fiuid with T^^^essed moving in 
the an inert spacetime gfj,u,inert- It is clear that the inert background is just the non-dynamical 
part of the metric, i.e., g^u,inert = (jf^u- We would like to claim that = T'^'^. This then 



^^We alert the reader again to footnote 18 - the expressions in (4.1) should be dealt with care as the l.h.s. 
and r.h.s contain contributions from quantities defined with respect to different metric structures. 
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would be a complete physical picture of how the dynamics on a Dirichlet hypersurface in the 
bulk emerges directly from the boundary description. 



4.1 Conservation equations at the boundary and on the Dirichlet surface 



Let us now implement the dressing picture heuristically described above at the level of the 
equations to derive the conservation equations on the Dirichlet surface from those on the 
boundary. Given an arbitrary energy-momentum tensor 



(4.2) 



with TT^'^ capturing the dissipative terms involving at least one gradient of the velocity field 
or thermodynamic state variables, we have 



V^T'^" = yV,e + {e + p)V,m1 + {e + p)a^ + P^^V.p + V, vr 



'"snd 



(4.3) 



where we have introduced 



■"snd 



dp _ dT 

de de ' 



and used the Euler relation e + p = sT , with s being the entropy density of the fluid. Since 
the part proportional to and the part transverse to should separately vanish, we get 



(4.4) 



Similarly, V^T'^'^ = is equivalent to the equation 



Using the relations 
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which are derived in §4.2 we can write 
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so that (4.5) becomes 



1 + 



d In a 
d\n Tj 



a 



^-^ ^ d In t) 



(4.5) 
(4.6) 



(4.7) 



(4i 



24 



For the equation (4.8) to describe the same dynamical system as the equation (4.4), it is 
necessary and sufficient that 
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1 + 



In T 



Ts 



P^jTVxI^/ -Cl^d^^ TT^^VaU^ 



Ts 



(4.9) 



We can show that this is indeed true at the first derivative level by using the conversion 
formulae (rewriting (3.27)) 




(4.10) 



d\na d 2 ^^ n ^ n 
dm I 2 a 

1 + ^{a^ - 1)^ ttu, 9u^ = 9u^ 

Csnd = U + 2*^" 

Hence, till this order, we have proved that T^^ is indeed the dressed energy-momentum 
tensor that we were looking for. 

It should be instructive to extend this analysis to higher orders in derivatives. In partic- 
ular, it would be interesting to pin down the specific property of the Dirichlet constitutive 
relation which leads to the fact that the dressed viscosity fj is same as the bare value t] 
and furthermore understand why the hypersurface fiuid has no bulk viscosity. This would 
complement the analysis of [34] who demonstrated the absence of corrections to the shear 
viscosity by considering a fiow equation in the linearized regime between the boundary and 
the horizon. 



4.2 Causality and relativistic fluids on the Dirichlet hypersurface 

Having established a clear connection between the dynamics of the dressed fiuid on the 
Dirichlet surface and that of fiuid on a 'dynamical' boundary metric, we now turn to ex- 
amining the properties of the fiuid motion. It seems a priori that all is well in the long 
wavelength regime with regards to the issues raised at the beginning of section §3 viz. the 
issue of locality and causality of the Dirichlet problem in AdS^+i. However, this is prob- 
ably a bit too quick; while it is true that we have local dynamical equations given by the 
conservation of the hypersurface stress tensor (3.24), we have not established firmly that 
these equations arise from a sensible thermodynamic system. We now proceed to address 
this issue. 
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The energy momentum tensor of the dressed fluid on the hypersurface is characterized 
by an energy density e and a pressure p which are given in (3.19). In particular, the pressure 
of the fluid is 



P 



[1 + 



a 



e which is useful in what follows. 



We also note that e 

Using the thermodynamic relations 

de dT dp 
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and e + f) = Ts 



e + p T e + p 
we get the entropy density and the temperature of this fluid as 
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(4.13) 



as quoted above in (4.6). The first of these relations follows from the fact that the entropy of 
the fluids on the asymptotic boundary as well as on the Dirichlet surface are given in terms 
of the area of the horizon which is unchanged by the solution. To determine the temperature 
on the hypersurface one has to account for the fact that the surface is in the interior of the 
spacetime. In the planar Schwarzschild-AdS^+i solution we get deviation from the Hawking 
temperature (which is temperature in the CFT) via a red-shift factor a. Conversely, given 
the above relations (4.13), the expressions for i and p can be deduced using de = Tds and 
dp = sdT . 

The speed of sound mode in this system is given by 
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de d-l 

This exceeds the speed of light as measured by g^i,, i.e., we get superluminal sound propa- 
gation, for a > ctsnd where ctsnd = ^3 — |. This corresponds to 
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^ 1 + - \n{3/2) + 0{d-^) 

One can intuitively understand this result from the viewpoint of the boundary fluid. As 
we noted earlier the boundary fluid is subject to a gravitational potential well. Should one 
locally increase the strength of this well then the fluid would get sufficiently accelerated, 
perhaps leading to a pathology. This is manifest in the picture of the dressed fluid moving 
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on an inert background achieved by translating over to the Dirichlet surface. In particular, 
this gets reflected in the fact that the effective pressure p felt by the dressed fluid increases 
relative to its energy density e thus driving the dressed fluid into a regime where the dominant 
energy condition is violated. Such violations of the dominant energy condition are known to 
be susceptible to superluminal sound modes^^ as we observed above. 

How much should one be worried by this apparent acausal behavior where the dressed 
sound mode travels superluminally with respect to the inert part of the boundary metric? 
After all, as is easily verified the mode with dispersion u ~ Csnd k propagates within the local 
light-cone of the 'dynamical' boundary metric g^^. This is achieved by the phenomenon we 
had already alluded to towards the end of the section §3.2.2: as we move our Dirichlet 
surface into the AdS, the Dirichlet constitutive relation ensures that the boundary light- 
cone opens up (see (3.28)) thus ensuring that the dressed sound mode is not superluminal 
when measured with respect to g^u- Of course, pending a detailed analysis of the initial 
value problem posed by this dynamical system (and other possible global issues), one cannot 
assert that the boundary physics is sensible from above observations alone. 

We would however like to suggest that viewing the hypersurface fluid as an autonomous 
dynamical system, a superluminal sound mode probably indicates a pathology. As described 
in [36] one should anticipate that the corresponding initial value problem^'' for the hypersur- 
face fluid might be ill-posed. Is this the way the long wavelength problem is telling us that 
the dual of the generic bulk Dirichlet problem in the CFT is ill-posed? Is it possible that the 
bulk Dirichlet problem in gravity is pathological the moment is finite and the fact that 
the effective dynamics of the fluid of the CFT remains sensible up to a critical radius rD^snd 
is just a long wavelength artifact? Clearly this issue deserves further investigation. 

We will now take an alternate approach which sidesteps these deep questions - given 
that our issue is with the superluminal sound mode on (for ru < ro^snd), is it possible 
to project this offending mode out of our dynamics hence avoiding the entire issue? We will 
now argue that fortunately the answer is yes - there is indeed a way to project out the sound 
mode, retaining sensible physics at least within the long wavelength regime. 

The way to do this is to move to the incompressible non-relativistic regime of fluid/gravity 
correspondence first studied in [13]. We will now implement their construction in our setting 
allowing us to obtain sensible dynamics for the Dirichlet problem. We will postpone some 
of the more general questions raised above to the discussion section §7. 

^■^For an early discussion see [35] where similar issues for fluids models driven to high pressure regimes are 
discussed. 

^^While the initial value problem for ideal, or even viscous fluids is ill-posed as the conservation equations 
are parabolic, we here want to drive home the point that the pathology we want to encounter happens for 
the sounds modes that are usually non-problematic. 
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5 The Dirichlet problem for gravity with non-relativistic fluids 

The discussion so far has concentrated on mapping relativistic fluid dynamics on a curved 
background S/j to a corresponding problem on the boundary where we have a relativistic fluid 
on a 'dynamical metric'. Moreover in the previous section, we cited the possible problems 
with the sound mode as a motivation to project it out consistently so as to get a clearly 
sensible dynamical system with no possible issues with the initial value problem etc.. Our 
goal now is to describe how this can be done consistently within our gradient expansion 
inspired by the non-relativistic incompressible scaling limit of [13, 23]. As we will see this 
limit has the added advantage that it naturally allows us to make contact with the metric 
derived in [22]. 

The idea as explained beautifully in [13] is the following: every relativistic fluid has a 
scaling limit where we freeze out the propagating sound mode, which drives the fluid into a 
non-relativistic regime, while simultaneously making it incompressible. This BMW scaling 
can essentially be derived by the requirement that one retains the non-linearities of the 
conservation equation (at least at first order). The resulting conservation equation is the 
classic incompressible non-relativistic Navier-Stokes equations. Using the fluid/gravity map 
[13] constructed a gravitational dual of this system. 

The BMW scaling involves two ingredients. Firstly, the velocities and the temperatures of 
the fluid are taken to be slowly varying functions of a specific kind, with spatial and temporal 
gradients having different scaling dimensions (heuristically dt ~ d1). Secondly, the amplitude 
of the spatial velocity and the temperature fluctuation (about some constant equilibrium 
value) are also taken to be small and are of the same order as dx and dt respectively. It is 
convenient to introduce a large parameter K (which is inverse of the small parameter e in 
[13]) in terms of which the above statements can be written as dt ~ dx ~ etc. We 
will denote the corresponding parameter for the hypersurface fluid by the hatted 

Under the large K limit it is possible to show that the relativistic conservation equations 
map straightforwardly into the incompressible Navier-Stokes equations. We review this 
scaling in Appendix D for the convenience of the reader and proceed in the main text to 
directly implement an analogous scaling on the hypersurface using a large parameter K. 

5.1 Non-relativistic fluids on the Dirichlet hypersurface 

To keep the computation sufficiently general we will take the metric on the Dirichlet surface to 
have non- vanishing curvature. Furthermore, it is useful as in [13] to allow for the background 
metric to be decomposed to slowly varying parts of different orders so as to recover non- 
relativistic fluids which are forced on E/). One reason for doing so is that we are going to 
obtain a boundary metric, via the map (fD described in §3.2.2, which naturally contains such 
terms. Hence it pays to be more general to see the mixing of various contributions at the 
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p) = -de + ^Jf (x) dx' dx^ (5.2) 



boundary. 

To start off let us consider on the hypersurface a metric g^^, of the form: 

9^.. = 9^} + V (5-1) 

with 

9 

where gf^^ are slowly varying functions of x* and with h^y are the metric perturbations which 
we take it as 

dx^ dx" = 2 k* dt dx' + (h*^ dt^ + h*j dx' dx^^ (5.3) 

To keep things simple it turns out to be useful to work with the background spatial metric 
gllj^ being Ricci fiat i.e., R^'^^ = 0. This turns out to simplify the analysis considerably for a 
host of terms dependent on the curvature of g^*^^ drop out - a more comprehensive analysis 
for general backgrounds is presented in Appendix E. We indicate the various corrections that 
arise from the curvature terms at appropriate stages in the main text. 

All the functions which have a * subscript or superscript (which we freely interchange to 
keep formulae clear) are of a specific functional form with anisotropic scaling of their spatial 
and temporal gradients. 

i(t,xO : M'^-i'i h^M , such that {aj,(t, a;^, vf^l(t, x^} ~ {©(«-'), 0(^-1)} (5.4) 

where K is a counting parameter introduced to implement the BMW scaling (on the boundary 
= eBMW as discussed in Appendix D). 

Following [13], we parameterize the velocity field as 

= (^1, Sf^-^ t):) (5.5) 

where the function is determined by requiring that g^yU^u'^ = —1. This gives the full 
velocity field in a large K expansion as 

u 



l + ^ihlt + l k* vi + gfj vi vt) + 0(K-^) 
= v: + ^ (h;, + 2 k* vi + ^ vi vl) vl + 0(K-^) 
= -1 - 2 {-Kt + 9$ vi t) + 0{y<'') 



= iv* + k*) + 



K,vl + \{h*u + 2k*vi+ ^ vi t) (v* + K 
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Note that most of these expressions are readily obtained by just 'hatting' the formulae in Appendix D 
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and the velocity gradients are given by 

e = 0(K-4) 

A^dx^ = a^dx^" = \dtv* + viwfv* - ff\ dx' + 0(K"^) 
a^^dx^'dx'' = k'^V^^v*^dx'dx^ -2k'HiV^^v*^dx'dt + 0{^<~^) (5.7) 

where V^"'' is the covariant derivative compatible with g^^^x"^) and we have freely raised and 
lowered the spatial indices with g^j^ for brevity. Further, /* is a forcing function determined 
as a functional of /i^,^ data 

f, = ^dih;,-d;k* + q%vi. (5.8) 

and q*j = vf^k* — vf^k* . In deriving these expressions we have used the fact that to 
leading order in the K — )■ oo expansion, the velocity field vl is divergenceless (see below). 

We take the scaling in b to be of the form 

b = bo + ^-^6K. (5.9) 

Using (3.20) 

a = ao + ^-' (Uo{l-&l)^-^)+0{^-') , ao^^=^== (5.10) 

"0 J VfioorD) 

we can evaluate the non-relativistic pressure per mass density and kinematic viscosity: 

6p 



£o+Po 



a// {bo td) 
lemat 

(i + ^ ("?«) - 1)) V 



e„ = ^ = ^ + o(«-) (5.11) 

with po = £o + Po = leTrc""! 6^ playing the role of the non-relativistic mass density. 

Given these data^^ we can show that the conservation equations (3.24) reduce to the 
incompressible Navier-Stokes equations:^'' 

Vf^p. + dtv* + vi Vfv* - 2 z>o V(°)^' {Vfv*^ = h (5.12) 

We now proceed to derive the expressions entering into the bulk metric and the map (pn 
given in §3.2.2. 

^''Note that BMW limit does not involve any scaling of the metric on E^i which is taken to be a fixed 
Lorentzian structure. Only the hydrodynamic fields (velocity, pressure etc.) are rescaled and their gradients 
are constrained to scale in a specific way. Thus the dressed fluid while non-relativistic still lives on a 
Lorentzian geometry. Further in this limit, the fluid dynamics has certain Galilean symmetries which are 
actually enlarged into a enhanced symmetry algebra [37, f3, 38]. 

^^If the spatial geometry on which the non-relativistic fluid moves is not Ricci-flat then there is an addi- 
tional term in Naiver-Stokes equations due to the background curvature, see (D.14). 
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5.2 Bulk metric in terms of Dirichlet data 

Armed with the results from §5.1 we can proceed to use those of §3.3 to construct the 
bulk metric corresponding to the non-relativistic fluid on the Dirichlet hypersurface Hd- In 
principle Einstein equations need to be solved in the new gradient expansion to obtain the 
non-relativistic solutions. As argued by the authors of [13], this can be done via an algorithm 
very similar to the algorithm used to flnd the metric dual of the relativistic fluid. The main 
difference is the anisotropic scaling of space with respect to time and the fact that the bulk 
metric is no-more ultra-local in space but is still ultra-local in time. We can again proceed 
from the space of non-relativistic solutions with asymptotic boundary conditions that we 
present in Appendix D and reparametrize it in terms of Dirichlet data. 

But we will take here instead an easier route and directly derive it from the the bulk 
relativistic metric (3.29) written in terms of Dirichlet data. We should be careful though 
- given the difference in derivative counting between the relativistic scaling and the non- 
relativistic scaling, it is in principle possible that a higher order term according to the 
relativistic counting contributes at a lower order according to the non-relativistic counting. 
In order to obtain the metric accurate to the order where the Navier-Stokes equations can be 
seen, 0(K^^), we need to have the certain terms in the relativistic metric accurate to third 
order in gradients. 

If we however restrict to the case of Ricci flat spatial metric on the hypersurface S^i, 
then we can obtain the non-relativistic metric from the second order relativistic fluid/gravity 
metric obtained in [12]. There are three terms we need to account for which give rise to non- 
relativistic contributions proportional to vf^V^^^^v* = V^q-j-O* and Vf''qii. These involve 
new radial functions which we collect below after presenting the bulk metrics highlighting 
the terms that were missed in the original analyses. General expressions including spatial 
curvatures can be found in Appendix E. 

Since in the scaling limit K ^ 1 one has from (5.7) that the = to leading order 
things simplify considerably. Using the formulae in the last subsection and including the 
terms from the second order metric^^ (highlighted) it is easy to show that the bulk metric 
(3.29) becomes^o 

^^This point was initially missed by both our analysis and that of [13]. It was originally thought that it 
would be sufficient to obtain the non-relativistic metric from just the first order relativistic metric. This 
unfortunately is not true and as a result the non-relativistic metrics quoted in vl of this paper and in [13] 
only solve Einstein's equations to first order in the non-relativistic gradient expansion. The correct form of 
the general expressions are now collected in Appendix E. 

^^Note that the highlighted terms involve Weyl invariant Ricci 3l^i, and Schouten curvature tensors 
which are defined in Appendix D. 

■^''Notc that we have retained certain terms at ©(H^"^) which are actually not necessary to solve the equation 
of motion (3.1) at this order (eg., the velocity cubed term). This is to facilitate ease of comparison of our 
results when we undertake the near-horizon analysis in §6, with those in the existing literature. 
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ds = —2a dx^dr + 
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The function L{x) which enters into the above formulae is given as 
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(5.16) 



while the other functions Mi and M2 can be found in (E.7). 

In the final results we have highlighted the terms appearing in ds\ that were missed 
in the first version of the paper and are necessary in order to solve Einstein's equations 
(3.1) to 0(K^^). Having realized the necessity of these terms it a-posteriori became clear 
that the BMW scaling metric introduced in [13] to describe the bulk dual of boundary non- 
relativistic fluids also receives corrections. The corrected form of this metric including effects 
of boundary spatial curvature is now presented in Appendix D. The general expressions for 
the Dirichlet non-relativistic fluid are collected in Appendix E as they involve many more 
terms compared to what was originally reported in vl of this paper. 

As required this metric reduces to the metric derived by the BMW scaling on the bound- 
ary [13] that we review in Appendix D in our conventions; we simply set r^) = 00, do = 1 in 
the above metric. As remarked above, it compares with the metric given in [13] up to the 
new terms mentioned above. 

Inspired by the gravity duals of fluids on cut-off hypersurfaces in flat space [16, 19], in 
[22] the result for the bulk dual to a non-relativistic fluid on a Dirichlet hypersurface in AdS 
has been recently derived. The results there are presented for a fluid on a flat background 
and are entirely contained within our framework. To facilitate ease of comparison with their 
results we now specify our computation to the case where the metric on Ej) is flat and 
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furthermore switch off the forcing. Taking g^^,, = r]^„ which amounts to setting glj — 
and k* = h^^ = h*j = in the above one can show that the bulk metric reduces to th 
(5.13) with 
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which agrees with that derived in [22] once one accounts for some differences in convention. 
In particular, one has to rescale the time coordinate to absorb the factor of do, in addition 



to redefining dt — )■ dg dt along with v* — )■ ^ v* and vl — )■ do-O*. The inhomogeneity in the 
scaling of the spatial velocities results from the fact that we define our hypersurface metric to 
be the induced metric rescaled by a factor of r^,^, while [22] works with the induced metric. 
Also, our derivation here allows us to go to 0(K~^) which is necessary in order to see the 
dynamical Navier-Stokes equation on the hypersurface S^i. 



5.3 The boundary data for the non-relativistic Dirichlet fluid 

Apart from the construction of the bulk dual to the non-relativistic fluid on the Dirichlet 
surface, we would also like to know what the corresponding physics on the boundary is. For 
instance we see that the boundary velocity field can be read off from the terms with dr in 
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(5.13) and is given as 
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while the boundary metric can be read off from the large r behavior and is given to be 
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with I/i(oo) denoting the asymptotic value of Liibr). 



Note that even in the non-relativistic limit the boundary 'dynamical' metric gfj,u^ light- 
cone is being opened up relative to that of the hypersurface metric. As before this is a 
consequence of the red-shift effect. We are normalizing here the hypersurface metric to have 
flat Minkowski metric to leading order and this causes a rescaling by an amount dg on the 
boundary. As long as do is finite this is not an issue for we are just encountering an overall 
rescaling of the boundary time. 

What we have derived here is the AdS analog of the membrane paradigm connection 
recently proposed in [16, 19]. Recall that the construction described in these papers proceeds 
by looking at an asymptotically flat geometry with Dirichlet boundary conditions at some 
timelike hypersurface (the analog of our S^)) and one solves vacuum Einstein's equations. It 
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was shown that the hypersurface dynamics is constrained to obey the incompressible Navier- 
Stokes equations, just as what we have shown above. However, the solutions described in this 
section solve Einstein's equations with a negative cosmological constant and we furthermore 
have argued that the Dirichlet dynamics is obtained by suitably dressing up of a CFT fluid 
by allowing it to propagate on a 'dynamical' background metric. In our context it is clear 
that the interpretation of the physics is less in terms of an RG flow, and more along the 
lines of the medium dependent 'dressing up' of the boundary fluid dynamics in contrast to 
the Wilsonian RG perspective put forth as an interpretation of the membrane paradigm 
originally in [15]. 

There is however one regime where our analysis should be able to make some contact with 
the discussion in [16, 19]; this is the near horizon regime where one expects to encounter a 
local Rindler geometry, which is the starting point for analyzing the Dirichlet problem in fiat 
space. In the next section we show how one can embed this construction using the solutions 
we have described, enabling one thus to explore the AdS version of the membrane paradigm. 

6 The near horizon Dirichlet problem 

So far in our discussion the Dirichlet hypersurface Hd has been located at some radial 
position ro that is finite. We now want to investigate what happens as we push this surface 
closer towards the horizon, i.e, S^i — t- "K^ via r^j — t- b^^. To understand the resulting physics 
we first have to realize that we are doing something strange: the horizon is a null surface 
and has therefore a degenerate metric. S/j on the other hand is constrained to be a timelike 
surface with a non-degenerate metric. So it is clear that demanding a well-behaved metric is 
going to result in an infinite scaling by the red-shift factor a; the main question is whether 
one can implement the scaling while retaining interesting dynamics. We will now proceed to 
show that there exists a scaling of parameters such that the near-horizon geometry makes 
sense. Furthermore, we argue that this allows us to embed the fiat space constructions of 
[16] into our AdS set-up. 

We outline the construction in a couple of stages to guide intuition: firstly in §6.1 we will 
examine the conservation equation on the hypersurface and from there infer the scaling of 
parameters. This is the only sensible thing to do for us, since the entire dynamical system 
of the boundary CFT fluid has been converted into that of the hypersurface fluid living 
on an inert background. We then examine the consequences of the scalings we derive in 
§6.2 focussing on the region close to the horizon; this amounts to blowing up the region 
of spacetime between 'K'^ and S/j. This blown up region bears close resemblance to the 
solution of the vacuum Einstein's equations dual to the incompressible Navier-Stokes system 
discussed in [16]. Indeed one should anticipate this based on the usual intuition that near 
any non-degenerate horizon one encounters a patch of the Rindler geometry. However, we 
will also encounter differences owing to the fact at the end of the day we are solving (3.1) 
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with a non- vanishing cosmological constant. 

There is also the further question of what the geometry between the Dirichlet surface 
and the asymptotic AdS region looks like and moreover what is the boundary physics of 
our scalings? We argue in §6.3 that the near horizon scaling regime renders the bulk metric 
viewed as a Lorentzian metric on a spacetime manifold nonsensical. However, it turns out 
that there is a nice language to describe the geometry in terms Newton-Cartan structures 
and we show that the bulk co-metric (which in usual Lorentzian geometry is the inverse 
metric) is well behaved as is the boundary co-metric. The result is quite satisfying from the 
physical perspective: the near horizon limit demands a drastic modification of the boundary 
metric, which forces one into a non-relativistic or Galilean regime. Consequently, rather 
than describing geometry in terms of Lorentzian structure, we are forced to use the less 
familiar but equally effective geometrization of the idea of a Galilean spacetime, in terms 
of Newton-Cartan geometry (see [39] for a nice account of this subject and [40] for a more 
recent review). 

6.1 Scaling of the Dirichlet dynamics in the near horizon region 

To understand the behavior of the fluid on the Dirichlet surface as we push S £, closer to the 
horizon, we first look at the conservation equation in this limit. It turns out that demanding 
non-trivial dynamics on the Dirichlet surface forces one into a scaling regime of the fluid, 
effectively making it non-relativistic. However, the scaling we encounter is not quite the 
BMW scaling [13] discussed earlier in §5 but a slightly modified version of the same. 

6.1.1 A new scaling regime 

Ignoring for the moment the fact for r£, < r^i ^nd we are supposed to be projecting out 
the sound mode to ensure subluminal propagation on the Dirichlet hypersurface, let us 
write the relativistic conservation equations and examine them as we zoom in towards the 
horizon. Consider then the conservation equation on the Dirichlet hypersurface (3.24); we 
have analyzed this from a generic viewpoint in §4.1 but for now we will focus on the truncated 
equations to second order by setting li^^ = —27] a^'^ (cf., (4.2)). Projecting these conservation 
equations parallel and transverse to the velocity we get (using rj as given in (3.16)) 

^d-l)u^^ + e + V.a^'^ = 

ad 

- + i ("^ - 1)) ^ - - ° 

We are going to try to analyze these equations in the limit when d — )■ oo. 

From (6.1) it is clear that if we insist on leaving the hypersurface data independent of 
a then we find that we have contributions at different orders that need to be independently 
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cancelled. The most constraining equation is the O(d^) term from the transverse equation 
which demands that the spatial gradients of b vanish. Then at 0(1) we have to kill the 
acceleration and have a non-trivial equation from the longitudinal part. At order we 
would need to kill all terms that show up with the shear viscosity. The upshot is that we 
are left with vacuous dynamics on the Dirichlet hypersurface should b and u'^ be 0(1) as 
a — )• oo. 

While this sounds a bit strange, a moments pause reveals that this is indeed what one 
should expect on physical grounds. The horizon of a black hole is a null surface (it is 
generated by null generators) and in the process of moving the Dirichlet hypersurface to the 
horizon, we are effectively doing an infinite rescaling (hence a oo) to bring the Dirichlet 
metric Qf^^ to be timelike and non-degenerate. Before we do such a rescaling however, we are 
in the ultra-relativistic regime as far as the horizon goes - in such a regime it is natural to 
expect that there is no dynamics, the fluid streams along the null generators and is effectively 
frozen into a stationary flow. 

It is then clear that in order to obtain non-trivial dynamics on the Dirichlet surface one 
has to scale the fields b and ii'^ in some fashion. The crucial question is whether there is 
any scaling that retains non-trivial dynamics; operationally demanding that we obtain an 
'interesting' non-linear equation. Consider the following scaling:'^^ 

b = ^b, + ^6b^, (l + 0(x"^), x^^l) , d = xd. (l + 0(x"^)) (6.2) 

X 

where the functions with subscript * have specific functional form with anisotropic scaling 
of their spatial and temporal gradients. Specifically, 

l(t, x') : M'^-i'i ^ M , such that {9tl(t, x'), d^t, x')} ~ {0(x-'), O(x-i)} (6.3) 

where we assign gradient weight 1 to the spatial derivatives and 2 to temporal derivatives. 
This is inspired of course by the non-relativistic BMW scaling discussed in §5.1 and as 
discussed there the 0(x^^) part of the velocity field is fixed by normalization.^^ 

However, there is a crucial difference from the BMW scaling; the leading term in the 
expansion of 6, is growing with x, which seems to be an issue. Nevertheless, it is easy to 
check that under this scaling (which admittedly is obtained by demanding a sensible 0(1) 
equation from the longitudinal part) one finds that the equations can be reduced to: 

7^2 I. 

d^vi = , -,^d,5b, + dtvt + vid.vt - -^Vht = (6.4) 
2o, a, a 

where we have specified to a flat Minkowski background metric for speciflcity [g^u = t]^^) 
and to make the equations more familiar. These are of course, the unforced (by background) 

■^^This scaling can be compounded with the scahng symmetry of Navier-Stokes equations [13] to change 
the exponents of but we refrain from doing so for simpHcity (see end of Appendix D). 

''^Roughly speaking i< ~ H of §5.1; this is the overall parameter that will organize for us the hierarchy 
necessary in the near horizon limit. 
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Navier-Stokes equations which we have earher described via the BMW scahng in §5.1. The 
incompressibihty condition as in that case arises from the longitudinal conservation equation 
at 0(x~^d) while the Navier-Stokes equation itself appears at (D{>i~^a). A similar scaling 
can be done for the forced Navier-Stokes system, we simply use the BMW scaling for the 
fluctuating part of the hypersurface metric. 



6.1.2 Deconstructing the near horizon scaling 

How do we reconcile this scaling derived here for the relativistic fluid with that derived by 
BMW [13]? Firstly, let's note that the rationale of our scaling up the background value of 
6 ~ is that shear term survives the scaling. From the second equation of (5.11) we 
clearly see that this is required in order to retain the shear term in the limit. On the other 
hand the non-relativistic pressure in the BMW limit is proportional to ^ and in the near 
horizon limit both bo ~ x6, and aQ k diverge. The extra scaling down of Sb^: l^b^, 
(over and above the scaling by is necessary to offset this divergence and ensure that the 
pressure gradient term contributes at the same order as the convective derivative dt + v^di 
and the shear. The effective pressure that enters into the equation (6.4) is 

X p^ = --xa,— — — (6.5) 
2 xo, 

using the scaling of a in (6.2) so that the first term of the Navier-Stokes equation in (6.4) 
is essentially dip^,. So the final equations of motion on the near horizon region for the 
hypersurface dynamics are simply: 

div\ = , dip. + dtvt + vidjv'l - z>.V'< = (6.6) 

with kinematic viscosity 

u. = 

a, d 



^. = A (6-7) 



We have here glossed over the fact that since we scale b k we potentially have a 
problem. The limit seems to suggest that we are taking the zero temperature limit of the 
black hole geometry since the Hawking temperature (which is seen on the boundary) is 
T ~ b~^ (4.13). This naively sounds like we are outside the long wavelength regime and as 
a result should not be using the fluid/gravity map to describe the Dirichlet problem. 

A different way to say this from a dynamical equation of motion perspective is that the 
scalings (6.2) were derived by demanding that the equations (6.1) remained non-trivial, which 
operationally means that different terms appear to have inhomogeneous weights. Hence by 
suitable fiddling of amplitudes and derivatives we engineered that terms which originally 
scaled as a to various powers all contribute homogeneously. However, we have not analyzed 
the higher order contributions to the equations of motion. Is it possible that relativistic two 
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derivative terms in the stress tensor, which show up as the corrections to viscous relativistic 
hydrodynamics equations of (6.1) show up at the same order? Note that this is not a problem 
for the case discussed in [13] for they engineered their scalings about a fixed background 
temperature, and there were no stray factors of K (equivalently k) in §5.1 (in front of Bq) to 
augment higher order contributions. 

For us to make the argument that the higher order terms are suppressed requires knowl- 
edge of the transport coefficients at second order on the Dirichlet hypersurface. This is in 
principle calculable by extending the Dirichlet fluid/gravity map of §3.2 to one higher order 
in gradients. However, we will now argue that there is an essential simplification which 
allows us to make the statement boldly without computation. The rationale is simply that 
the enhanced scaling of the zero mode part of b is compensated by the suppressing the cor- 
responding fiuctuation term 6b so as to ensure that is finite. Since the overall value of 
bo scales out in the non-relativistic regime (we will see this clearly from implementing the 
scahngs in the next subsection) it cannot affect the resulting equations. So we conjecture 
that accounting for higher order corrections as well, one will obtain (6.6) as the leading order 
equations in the x expansion (all further corrections are suppressed by higher powers of x). 

One physical way to motivate the correction is to first note that while the boundary 

temperature is being scaled to zero in the limit (6.2) the hypersurface temperature can 
indeed be maintained to be finite:'^'^ 

T. = -V 6.8 

So from the Dirichlet observer's point of view there is still scope for a non-relativistic scaling, 
and in fact this is how the Dirichlet observer would carry out the BMW regime. What is 
clear is that due to the extra red-shifting in translating to the boundary, the asymptotic 
observer is going to have trouble reconciling this near horizon scaling regime in his variables. 
We will return to this issue in §6.3. 

In summary, the non-relativistic incompressible Navier-Stokes equations determine the 
dynamics of the hypersurface fiuid as S^i approaches the horizon; we will now proceed to 
investigate what this means for the various metrics: first we look at the bulk metric first and 
then examine what is the effect on the boundary. 



6.2 The bulk metric between E/j and JC"*" 



From the discussion at the end of §6.1.2 it seems natural that we should try to ask the 
question as to whether we can satisfy the constraints of the Dirichlet problem in the re- 
gion between and the horizon (for the moment forgetting about the asymptotic region). 
Clearly we have T^d approaching the horizon, so we will be forced to consider a double- 
scaling regime where we zoom in close to the horizon and expand out the spacetime region 

■^■^We thank Sayantani Bhattacharyya for emphasizing this point to us. 
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in-between. Let us temporarily ignore the consequences of this on the equations of motion 
(3.1) and formally take the limit of the bulk metric given in (3.6), (3.29). 

To zoom into the region between the horizon and S^, realizing that r^i — )■ | in the limit, 



we take: 



from which it follows that 



6, X V k%, a. 



dr = . dp 

K-^ 0, a. 



Similarly we parametrize the position of our Dirichlet hypersurface via 

-^ = r^fi + -^V (6-10) 



6, X \ x^fe, d. 

Both py) and a, are a measure of how close we are to the horizon after one has zoomed into 
the near horizon region. To relate them, we subtitute into the expression for a and do a 
large x expansion. Identifying the leading term as d,, we get 

Pd = - — -, = t'. (6.11) 
a, d 

We see that p_D is same as the kinematic viscosity - hence, the po limit is identical to 
the inviscid limit in hydrodynamics. 

Implementing this change of variables and the scaling (6.2) in the bulk metric we obtain 
the near-horizon metric. In fact, the fastest way to derive the metric in (6.14) is to utilize 
the fact that the near horizon limit is essentially the non-relativistic limit together with a 
particular form of the pressure fluctuations; essentially we substitute do = xd,, bo = x6,, 
p^: = into (5.17) along with the identification K = x to ensure that we have the correct 
scalings. We first calculate the near horizon expansions of various functions appearing in the 
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metric 
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Using these, we get the metric as 
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(6.13) 



df 



2 4 a ai 



dx" dt 
dp dx^ 



~ [P^ - Pd + ^ Pd {pd - p)] djd^v* dx^ '^^ + ^ djdH'l dp dx^ 

(0) 



(6.14) 



We have restricted attention to the simplest setting where g^y = rj^u (hence g\j = 5ij above) 
and -0^ = Sijvlvi. In deriving the expressions as in §5.2 we have used the incompressibihty 
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condition, which appears as before as the leading order fluid equation of motion. The Navier- 
Stokes equations themselves can also be used to simplify the last term in (6.14): one can 
replace the convective derivative of the spatial velocity dtv* + vldjv* by —dip^, + 2 0, V^f *. 

This is in fact closely related (but not identical) to the metric derived in [16] and gener- 
alized further in [19] to higher orders (compare with for instance Eq (3.2) and (6.5) of the 
latter paper). We recall that these works consider solutions to vacuum Einstein's equations 
without a cosmological constant. In particular, [16] solves the Dirichlet problem in flat space 
by looking for small gradient fluctuations around a Rindler geometry, while we are dealing 
with solutions of (3.1) which has a non- vanishing negative cosmological constant. 

We see that to leading order we nevertheless should expect agreement to the analysis 
of [16, 19] - this is clear from the first two lines of (6.14). This is reflecting the universal 
Rindler nature of a non-degenerate horizon such as that of the planar Schwarzschild-AdS^+i 
solution we started with. Moreover, we can also physically understand how a solution of 
(3.1) reduces to that of vacuum Einstein's equations by noting that in the limit we consider 
here, the cosmological constant gets diluted away. The factor of k on the l.h.s, is essentially 
indicative of -RacIS RAds in our near horizon limit. This is to be expected; by zooming 

in onto the region between T^d and we are effectively blowing up a small sliver of the 
spacetime and are thus losing any information about the background curvature scale -RacIS- 
In some sense the near horizon limit is like the limit we take to decouple the asymptotically 
flat region from the AdS throat in the D3-brane geometry; what is unclear is whether some 
notion of decoupling exists in the present context. 

But starting at 0(x~^) on the l.h.s. of (6.14) we start to see deviations from the metric 
presented in [16, 19]. Even though we are zooming close to the horizon of a Schwarzschild- 
AdSrf+i black hole, starting at second order we should expect to see curvature contributions. 
The various terms at this and higher orders originate from how the Rindler geometry gets 
corrected as we step away from the strict limit. In particular, using the near-horizon scaling 
of variables we see that (3.1) reduces to 

RaB + 4 ^2 9aB = 0. (6.15) 

It follows that all terms which involve 4j in the denominator correspond to the corrections 
due to the AdS curvature. Setting such terms in (6.14) to zero i.e., taking a, — )■ oo leads us 
to the metrics derived in [16, 19]. 

We are here restricting attention to the metric to leading orders in the k expansion; upto 
0(x~^). In principle since the fluid/gravity metrics are known to higher orders one can carry 
out the construction to higher orders and we expect that we should be able to reproduce 
the results of [19] who have derived the expressions in the near horizon limit to 0(x~^) in 
our notation. Note that this simplification of the near horizon region construction happens 
because of the correlated amplitude and gradient scaling, as already emphasized in [13]. 
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6.3 Near horizon limit and Galilean degenerations 

Having seen that it is possible to take the near horizon scaling limit and retain both inter- 
esting dynamics on and further have a sensible geometry between Hd and CK"*", our next 
goal is to address whether we can make sense of the metric between the hypersurface and 
the boundary. There are already many indications that this is going to be problematic; the 
vanishing of the boundary temperature oc being the most prominent one, which suggests 
break-down of the gradient expansion. 

The bulk metric turns out not to make sense, but we shall see that there is an object that 
does - this is the inverse bulk metric which we shall refer to as the co- metric (see below). 
Similarly, the boundary co-metric also is well behaved and we will argue naturally provides a 
Newton-Cartan like structure on the boundary so that the boundary geometry degenerates 
from a Lorentzian manifold into a Galilean manifold. 

In the following, we will call the metric on the co-tangent bundle g'^^ as the co-metric. 
This is a more accurate terminology in the context of non-relativistic limit (the Galilean 
or Newton-Cartan limit) than the usually preferred 'inverse metric' since in this limit the 
co-metric degenerates (becomes non-invertible) and the usual metric ceases to exist. Hence, 
we will prefer as much as possible to work with the co-metric instead of the metric. 

6.3.1 Emergence of Galilean structure on the boundary 

Let us first examine what happens to the boundary metric data when we try to push the 
Dirichlet hypersurface towards the horizon. This corresponds to a tending to infinity. We 
will work with the relativistic expressions of §3 to maintain covariance; it is easy to then pass 
over to the near horizon scaling regime discussed in §6.1 and obtain explicit parameterization 
of the results. 

From the formulae for the boundary metric in terms of (4.1) we see that there are 
various terms which blow up. Moreover, in some cases the higher derivative terms overwhelm 
the lower derivatives suggesting a breakdown of derivative expansion, as already suspected 
from the vanishing of the boundary temperature. It is clear that this limit if it exists cannot 
be straightforward. 

As the Dirichlet hypersurface tries to approach the horizon, it first hits the hypersurface 
fD = fD,snd- From the boundary viewpoint, the interaction between the fluid with its 
gravitational potential packet drives the system to have superluminal sound propagation, 
Csnd exceeds the dressed speed of light as determined by g^'^. By this point the interaction 
between the boundary metric and the fluid has become so important that the bare velocities 

have no more physical significance. We should rather think in terms of the dressed fluid 
and see how we can make sense out of the approach to the horizon. 
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So, we will first rewrite the above formulae in terms of the dressed velocity m^. The 
dictionary for the (normalized) metric and the co-metric are 
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td [1 + 



We have already remarked on the necessity to project out the dressed sound mode and 
take an incompressible limit of the dressed fluid a la BMW as the Dirichlet hypersurface 
crosses ro = rD,snd in §5. The formulae above reinforce that intuition - we see for example 
that unless 6 is sent to zero, the first derivative terms in the boundary metric overwhelm 
the zeroth order answer thus leading to a breakdown of the derivative expansion. So, we 
will drop the 6 terms with the understanding that we are projecting into the incompressible 
sector of the dressed fluid. Even this does not seem to help the case for the metric, since 
it still diverges - keeping only leading order terms (and assuming none of the subsequent 
terms grow with d - we will postpone a more careful analysis for the future) we get 
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(6.17) 



We are thus led to a remarkable conclusion: the near horizon Dirichlet constitutive rela- 
tion is that the boundary co-metric degenerates along = i.e., g^^t^, = and the bound- 
ary metric has one divergent time-like eigenvalue along the same direction g^i, = —a^t^tv 
This is the signature that the boundary metric is becoming Newtonian/non-relativistic - 
crudely speaking this is analogous to the c — i- oo behavior of Minkowski metric/co- metric 

r]'"'' = diag(-^, 1, 1, . . .) and r^^,, = diag(-c^ 1, 1, . . .) 

The way to make this mathematically precise is to resort to what is called as Galilean 
manifold. 

A Galilean manifold is a manifold with a Newtonian time co- vector t^, and a degenerate 
co-metric g^'^ satisfying (/^^tj/ = 0. As in standard differential geometry, we can demand a 
connection that is compatible with this structure. This requires that there exists a torsionless 
connection defining a covariant derivative Vj^*" which is compatible with both the co-metric 
and the time co- vector, i.e., V^'^g^^ = and V^'^tu = 0. If in addition we have such a 
a torsionless covariant derivative compatible with the Galilean structure then we call 
V'^'" as the Newton-Cartan connection and the corresponding Galilean manifold is termed a 
Newton-Cartan manifold. See [39, 40] for further discussions and [38] for another AdS/CFT 
perspective on the Galilean structures. 
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Thus we have just argued that that there is a natural emergence of Gahlean structure 
in the boundary theory as we take our Dirichlet surface nearer and nearer to the horizon. 
In this case, the metric that we have to put the field theory in gets enormously simplified 
and we get just a CFT on a Galilean manifold with a Galilean co-metric having as the 
degenerate direction. This is a very precise way of stating what the membrane paradigm is 
from the boundary theory viewpoint - the membrane paradigm is a particular Galilean limit 
of the field theory. '^"^ 

Having shown that we have a Galilean structure on the boundary, we further ask: "Is 
there a Newton-Cartan structure" ? Consider the Galilean limit of the Christoffel connection 
which is by construction torsionless and compatible with the Galilean co- metric g^'^ . 
This could be a sensible Newton-Cartan structure if it annihilates the time co-vector; to see 
what V ^ty is we use 
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(6.19) 



This in fact determines a sensible answer in the infinite a limit , once we factor in our earlier 
observation that 9 is suppressed in the near horizon limit (it is 0(q;~^)) as are any other 
gradients of the velocity field. These scalings can be read off essentially from §5 as the near 
horizon limit is the BMW limit as far as the velocities are concerned. 

The leading contribution to (6.19) comes at 0(d;~^) from the vorticity and expansion. It 
is tempting to speculate that the degeneration of the co- metric along with (6.19) define a 
Newton-Cartan like limit and we will postpone a structural analysis of this limit for future 
work. In the next subsection, we will try to examine the bulk co-metric to see what this 
Galilean limit entails in the bulk. 

■^^Once we include the new corrections highlighted for e.g., in (5.17), we find that the boundary Dirichlet 
constitutive relations and the boundary gradient expansion are in tension. For instance we find terms which 
originate at 0(H~^) migrate down to 0(H~^). However, to the order we have looked at all such terms are 
related to enforcing Landau frame choice both on the hypersurface and the boundary, and hence to a gauge 
choice for the bulk metric. Our statements in this section assume that such pure gauge terms are irrelevant, 
but this issue deserves further investigation. 
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6.4 The bulk co-metric in the Near-Horizon limit 



In the preceding subsection, we argued for the emergence of Galilean structures in the Near- 
Horizon limit thus forcing us to formulate the boundary geometry in terms of a co-metric 
than a metric. Since the boundary metric is ill-behaved in this limit, it is clear that the bulk 
metric should be ill-behaved too. This poses a conundrum since it seems as if in taking the 
near-horizon limit we have destroyed the AdS-asymptopia and one might wonder whether 
there is any sense in talking about the bulk geometry far away from the horizon. We will 
in this subsection take the description of the boundary geometry via a co-metric as a clue 
and argue that the bulk geometry is also well-described by a bulk co- metric. Thus while the 
metric description might break down the co-metric description with its associated Galilean 
structures continue to describe the bulk and the boundary geometry. 

Now we present the co-metric as a function of the Dirichlet data. 
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As we anticipated, there are no problems evident in the infinite a limit (provided one takes 9 
to zero appropriately to preserve the validity of derivative expansion). Contrast this with the 
bulk metric (see (B.14)) whose large a limit seems dubious. This supports our contention 
that there is a completely sensible description of the bulk and the boundary geometry in 
terms of co-metrics everywhere when we take the near-horizon Dirichlet problem along with 
an incompressible limit on the dressed fluid. 

Our exercise of rewriting the near-horizon bulk metric in terms of the non-relativistic 
variables can be repeated in the case of co-metric and we can easily convince ourselves that 
this goes through without any new subtleties. In fact, most terms in the above expression 
are sub-dominant in the k expansion introduced in the previous subsections. In fact, at 
the zeroth order the spatial part of the co-metric is just that of vacuum AdS; with the 
temporal part contributing only at a higher order. The horizon structure encoded in {br)~'^ 
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is completely invisible until quite high orders in k expansion, which is to be expected since 
the boundary temperature is being scaled to zero. To be specific using the scaling (6.2) we 
find ^ 

g^^a^ ®dB = r^dr®dr^ -di®d,^ k-^ (2 dt ®s dr) + 0{k'^) . (6.22) 

Basically we seem to find that the co-metric description is mostly oblivious to the near- 
horizon geometry (which as we have seen is well-described in the metric description). Hence, 
we are naturally led to an effective description where there are two regions of the geometry: 
one well-described by a metric, and another by a co-metric. 

The interesting question is to see whether there is an overlap region where the two 
descriptions are valid and the metric is just the inverse of the co-metric. Having detained the 
reader for so long, we will leave the detailed answer of this question to future work. However, 
we would like to draw the attention of the reader to the following fact - the near-horizon 
limit of the Dirichlet problem naturally seems to lead to novel geometric structures closely 
associated with the Galilean limit. These Galilean structures evidently call for more detailed 
studies especially since they might hold valuable lessons for non-relativistic holography (as 
previously pointed out in [38]). 

7 Discussion 

The bulk Dirichlet problem in AdS, which we defined as the gravitational dynamics in a 
spacetime with negative cosmological constant, subject to Dirichlet boundary conditions on 
a preferred timelike hypersurface S^i, is interesting in the AdS/CFT context for several 
reasons. For one, it allows us to investigate questions in a cut-off AdS spacetime which 
might be dual to a field theory with a rigid UV cut-off via the usual AdS/CFT dictionary. 
Furthermore, it allows us to touch upon the ideas involving holographic renormalisation 
and the implementation of the Wilsonian RG ideas in the gravitational description. The 
main motivation behind our work was to ask, what is it that one is doing on the dual field 
theory living on the boundary of the AdS spacetime that ensures these Dirichlet boundary 
conditions on E^), and whether such a problem is always well-posed. 

For linear systems in AdS it is easy to see that there is a simple map between data on 
the hypersurface and that on the boundary. In particular, Dirichlet boundary conditions on 
S^) translate into non-local multi-trace deformations of the dual field theory. One can think 
of the boundary theory being deformed by some 'state dependent' sources. Such a boundary 
condition while seemingly bizarre from standard field theoretic perspective, has an effective 
description in terms of a much simpler and local source, which is the Dirichlet data on S^i. 

In this paper we have further argued that such bulk Dirichlet problems for full non-linear 
gravitational dynamics in AdS are also amenable to solution in a certain long wavelength 
regime. Using results from the fluid/gravity correspondence [8, 12] we have constructed 
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the bulk spacetime resulting from the specification of a hypersurface metric on S/). The 
solution is contingent upon the hypersurface dynamics, as determined by the conservation 
of the stress tensor induced on S^), be conserved. In the long wavelength regime this stress 
tensor can be written as that of a non-conformal fluid propagating on the fixed geometry 
of S^). Armed with such a solution we can examine the dual CFT to determine how one is 
achieving the Dirichlet dynamics. 

In the long wavelength regime the boundary theory is described by fluid dynamics with 
what we term as Dirichlet constitutive relations. In particular, while the stress tensor of 
the boundary fluid is of the conventional conformal fluid form, it lives on a background 
geometry whose metric is 'dynamical' in the following sense: the boundary metric depends 
on the dynamical fluid degrees of freedom. One should think of this in terms of a fluid being 
subject to a gravitational potential well which is furthermore carried along with the fluid. 
A useful analogy is to think of a charge carrier in a polarizable medium. The medium or 
more precisely in our case the 'dynamical' background exerts a force on the fluid. However, 
one can subsume the 'dynamical' background's effect on the fluid, and rewrite the dynamics 
as that of a 'dressed fluid' on a fixed background. This 'dressed fluid' is a collective effect 
and moreover from the geometry we know what its description should be - it is simply 
the non-conformal fluid on the Dirichlet surface. This can be independently verified by 
starting with the boundary fluid and the 'dynamical' boundary metric and showing that the 
boundary conservation equation can be rewritten as the hypersurface conservation equation, 
thus deriving the non-conformal fluid stress tensor on S/j. 

We also note that in the long wavelength regime because one is working order by order 
in boundary gradients, the boundary fluid is deformed locally - at any spacetime point 
the source g^i, depends only on the fluid degrees of freedom (velocity and temperature) at 
that point. This is what allows us to think of the fluid as carrying around with it a local 
gravitational cloud. From a formal point of view, the gravitational Dirichlet problem involves 
turning on local multi-trace deformations for the field theory. 

From the above discussion one is then tempted to say that the bulk geometry provides 
a way to repackage non-local deformations into a local perturbation at a lower radius or 
scale. This is highly suggestive of some sort of renormalisation of sources as one propagates 
boundary conditions into the bulk. Let us compare this picture with the recent discussion 
of the holographic Wilsonian RG flow idea of [6, 7]. In these works starting with a field 
theory with given sources for say just single trace operators in the planar theory, one derives 
an effective action containing not just renormalized single trace sources, but also multi- 
traces on some chosen cut-off surface. The flow equation governing the radial evolution of 
such sources was argued to arise by effectively integrating out a part of the bulk geometry 
between the boundary and the cut-off surface. It was also important for that discussion 
given the particularization to the Wilsonian perspective that one does not prejudice oneself 
with the boundary conditions in the interior of the spacetime below the cut-off (the infra-red 
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of the field theory). 

In the current context we see that by transferring the boundary conditions onto the 
Dirichlet surface we have shifted the burden of multi-traces onto the boundary - as de- 
scribed above there is a suitable set of non-local multi-trace deformations on the boundary 
which 'renormalizes' into a single trace source on the cut-off. Whilst this is not totally sur- 
prising for linear systems such as the scalar problem discussed in the text, it is satisfying 
that a similar statement can be made for non-linear gravitational dynamics by invoking the 
long wavelength gradient expansion. We should however note that our discussion explicitly 
assumes knowledge of the interior (infra-red) boundary conditions; in the fluid/gravity solu- 
tions constructed herein we have demanded that there be a regular future event horizon to 
single out a sensible solution. 

Having obtained the effective dynamics on the hypersurface, we learnt that the conser- 
vation equations result in a sound mode which travels outside the light-cone of the fixed 
hypersurface metric once one pushes the surface too far into the interior. Past this sonic 
threshold, it is probably not sensible to maintain a relativistic fluid description on E^- We 
argued that we should pass over into a non-relativistic regime, by a suitable scaling of fluid 
variables (a la BMW [13]). It is interesting to ask whether the acausality manifesting itself 
in the sound mode of the hypersurface dynamics can be discerned (without invoking the 
'dressed' picture) from the boundary. More importantly, one wonders whether the gravita- 
tional Dirichlet problem suffers from such pathologies in general. "^^ Does one encounter any 
acausal behavior the moment the hypersurface is inside the bulk outside the long wavelength 
regime? What about their validity where say stringy effects are taken into account? Can we 
engineer such a boundary condition consistently using objects like orientifolds in string the- 
ory? What happens once we go beyond large N? Does the field theory with these non-local 
deformations make sense beyond large N for at least some subclass of these deformations? 
What is the field theory interpretation of other kinds of boundary conditions other than 
Dirichlet boundary conditions (say if we fix the mean extrinsic curvature as proposed in 
[43]) ? These are fascinating questions which deserve to be explored further. 

Finally, by examining the Dirichlet problem in the vicinity of the event horizon of the 
fluid/gravity solutions we made contact with the recent constructions of flat space gravita- 
tional duals to incompressible Navier-Stokes flows on a cut-off hypersurface [16, 19], deriving 
in effect the membrane paradigm from the boundary field theory. Focussing on a tiny sliver 
of the geometry between the horizon and S^i leads to the long wavelength solution around 
the Rindler horizon found in these works. This is reminiscent of the Penrose limit; while in 
the latter we focus on the geometry close to a null geodesic and blow it up, here we zoom 

■^'^We believe the answer to this question is most Ukely to be in the affirmative. If we consider deforming the 
boundary conditions at infinity by making the boundary metric an arbitrary local function of multi-traces 
of the stress tensor, we can think of the problem effectively as a mixed boundary condition a la [41], which 
has been argued to suffer from ghosts generically (see also [42]). 
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close to a null surface and blow up the spacetime in its vicinity (perhaps a better anal- 
ogy is the near horizon limit of black D3-branes without any statement about decoupling). 
The process of blowing up the near horizon region dilutes away the cosmological constant; 
hence rather than obtaining a solution to Einstein's equations with negative cosmological 
constant we end up with a geometry that solves vacuum Einstein's equations at leading and 
next-to-leading order. However, starting at higher orders we start to see the effect of the cos- 
mological constant, revealing the throat region between the near-horizon Rindler geometry 
and the asymptotic AdS spacetime. 

From our embedding of the construction of [16, 19] into the fluid/gravity correspondence 
[8, 12] we learn that the boundary fluid lives on a manifold with Galilean/Newton-Cartan 
like structure in the near horizon limit. This degeneration of the Lorentzian structure into 
an effective Galilean one is what enforces the incompressible Navier-Stokes scaling from the 
viewpoint of the CFT. One can simply say that the membrane paradigm is a particular 
Galilean limit of the field theory. This perspective also clarifies the universality of the 
membrane paradigm - since we are zooming down to the Rindler region in the vicinity of a 
non-extremal black hole horizon, we should expect the same geometry for any non-degenerate 
horizon. Hence the dual description should always be in terms of an incompressible Navier- 
Stokes fluid as long as we deal with systems carrying only conserved charges [13]. In the 
presence of other light degrees of freedom, as happens in systems with spontaneously broken 
symmetries, say the holographic superfluids discussed in [45, 46, 47, 48] we would have to 
project out all linearly dispersing modes to achieve the same. 

In the near horizon limit we find that both the bulk metric and the boundary metric 
degenerate into a Newton-Cartan time-metric, but the co-metric is spacelike and well be- 
haved. We speculate that the bulk spacetime in this limit should be described in two patches: 
the region between the horizon and the Dirichlet hypersurface enjoys a description in terms 
of a conventional metric while the region between and the boundary requires use of 
the co-metrics and Newton-Cartan structures. We conclude that the language of Galilean 
geometries is what is necessary to implement the membrane paradigm in the AdS/CFT 
correspondence. 

The idea of implementing a Galilean limit of AdS/CFT correspondence was proposed 
recently in [38] who were also motivated by considerations involving the incompressible 
Navier-Stokes equations, which turn out to enjoy an enhanced symmetry algebra [37]. It 
was argued there based on this algebra that appropriate Newton-Cartan manifold should 
correspond to an AdS2 slice of the bulk AdS^+i spacetime. Here in contrast we are retaining 

■^^It has been pointed out in [44] that there is an obstruction to finding solutions with Dirichlet boundary 
conditions in the near horizon Rindler like region, when the fluids live on compact curved spatial manifolds, 
such a a sphere (as would be relevant for the near horizon of asymptotically flat Schwarzschild black hole). In 
our construction the near-horizon limit described in §6 requires that we simultaneously scale the curvatures 
to zero, and so we are unable to see the origins of such obstructions. 
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the entire manifold, not just the two dimensional slice. However, we achieve this at the 
expense of discarding the Lorentzian metric structure, and work with a co-metric and a time 
CO- vector. 

Clearly the role of Galilean/Newton-Cartan structures in the context of AdS/CFT re- 
quires further investigation. The bulk geometries we construct here seem to provide an 
interesting interpolation between Lorentzian and Galilean structures somewhat reminiscent 
of the Galilean limit procedure described in [49, 50]. It is also tempting to contemplate the 
possibility that apart from the context discussed herein, these structures could also provide 
useful clues on how to work with the gravity duals of non-relativistic field theories. In the 
context of Schrodinger spacetimes which were proposed as duals to non-relativistic conformal 
field theories [51, 52] the Galilean limit is usually implemented as a DLCQ limit which can 
naturally be studied in the Newton-Cartan language. 

There are many directions in which the constructions we have described here can be 
generalized. For one it would be very useful to understand the gravitational Dirichlet problem 
outside of the long wavelength regime and to investigate its well-posedness. Further it would 
be useful to flesh out in precise detail the connections (if any) between the ideas around 
holographic RG flows and the Dirichlet problem and ask whether one can use the ideas 
developed here to implement the Wilsonian holographic flow of [6, 7] in the non-linear gravity 
context. 

Even within the long-wavelength regime there are tantalizing similarities with the black- 
fold approach [53, 54], where one considers hypersurfaces which have intrinsic as well as 
extrinsic dynamics. Freezing out the extrinsic dynamics should lead one to the gravitational 
Dirichlet problem (for instance in the analysis of the black string and membrane instabilities 
[55] the extrinsic dynamics was frozen and the blackfold equations are precisely those of fluid 
dynamics). While the blackfold analysis is generically well suited for co-dimension three of 
higher hypersurfaces, it appears that one could recover some of the results discussed herein 
by extrapolating the equations to a co-dimension one hypersurface.'^^ 

Of more immediate interest is to complete the derivation of the hypersurface dynamics 
from that on the boundary to higher orders in the gradient expansion. In particular, while we 
have shown that the ideal fluid equations on S/) arise from those on the boundary, moving to 
second order in gradients in the conservation equation will allow us to show at the non-linear 
level that the shear viscosity of the hypersurface fluid is the same as that on the boundary. 
This statement would establish the non-linear version of the non-renormalisation of r] (or 
equivalently rj/s since we know that the entropy density being associated with the horizon 
remains unchanged) complementing the earlier analysis of [34] who showed this in the regime 
of linearized hydrodynamics using a flow equation. One would hope to also understand by 
this study why the bulk viscosity term on the hypersurface vanishes despite the fluid being 
non-conformal and having a non-vanishing trace for the stress tensor. At the same time it 

'^''We thank Roberto Emparan for useful discussions on this issue. 
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would be interesting to understand how the higher order transport coefficients evolve from 
the boundary to S/j which should be possible to examine within our framework with a little 
bit of work. 

Another interesting avenue for exploration is to ask about the gravitational Dirichlet 
problem in the presence of degenerate horizons. All of the discussion in the present pa- 
per and indeed in the fluid/gravity correspondence has focussed on situations of thermal 
equilibrium at non-zero temperature and hence one naturally studies spacetimes with non- 
degenerate horizons. Degenerate horizons pose an intriguing challenge and it would be useful 
to understand how the boundary theory (which should be more than just a fluid dynamical 
system) reorganizes itself as we pass over to the hypersurface description. It would also be 
interesting to extend our considerations to stationary black holes where it would be useful 
to understand the Dirichlet dynamics when the hypersurface is inside the ergorsurface (as 
defined by the asymptotic observer). 

Finally, another interesting avenue for contemplation is whether the ideas discussed herein 
can be ported to the context of brane-worlds, where we have gravity induced on a cut-off 
surface in AdS. Recently [56] attempts to derive the brane gravitational equations by working 
on a hypersurface near the boundary of AdS and invoking the Fefferman-Graham expansion. 
It is clear from our discussion that no matter what boundary condition we choose on the 
hypersurface, it is most likely to involve non-local deformations on the dual boundary field 
theory. What this implies for induced gravity scenarios is an issue that deserves further 
study. 
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A Notation 



We work in the ( — 1-+ . . .) signature. The dimensions of the spacetime in which the conformal 
fluid hves is denoted by d . We usually assume d > 2 unless otherwise specified. In the context 
of AdS/CFT, the dual AdS^+i space has d + 1 spacetime dimensions. The lower-case Greek 
indices ^,1^ = 0,1, . . . ,d — l are used as boundary space-time indices, whereas the upper-case 
Latin indices A,B = 0,1, ... ,d are used as the bulk indices. The lower-case Latin indices 
i = 1, . . . ,d — 1 index the different spatial directons at the boundary. 

Throughout this paper, we take the extra holographic co-ordinate to be r with the bound- 
ary of the bulk spacetime at r — )■ 00. 

Among the objects carrying Greek indices the hatted variables belong to a hypersurface 
r = r^) whereas the unhatted objects naturally belong to the boundary r = 00. Further we 
use * to mark non-relativistic objects. 

We use round brackets to denote symmetrisation and square brackets to denote anti- 
symmetrisation. For example, = ^ {B^,, + B^^) and = \ {B^^, — B^^). For tensor 
products we denote symmetrisation with an explicit subscript as X^^F = '^{X®Y + Y®X). 

A.l Quick reference table 

We have included a table with other useful parameters used in the text. In the table A.l, 
the relevant equations are denoted by their respective equation numbers appearing inside 
parentheses. 

A. 2 Notation in the Bulk 

We denote the bulk metric by Sab- The inverse metric in the bulk (we will call this the 
co-metric - since it is the metric on the cotangent bundle) is denoted by S^^- 

We take the bulk AdS radius to be unity which is equivalent to setting the bulk cosmo- 
logical constant to be A^+i = — We denote the bulk Newton constant by Gd+i- For 

the ease of reference, we now give the value of Gd+i for some of the well-known CFTs with 
gravity duals (see [1, 57, 58] for further details): 

1. The d=4, !N=4 Super Yang-Mills theory on A^^ D3-Branes with a gauge group SU(Nc) 
and a 't Hooft coupling A = Qym^c is believed to be dual to IIB string theory on 
AdSsxS^^i with G5 = Tr/i2N^) and a' = {AttX)-^/\ 

2. A d=3, 7^=6 Superconformal Chern-Simons theory on Nc M2-Branes with a gauge 
group U(Nc)fcX U(Nc)_fc (where the subscripts denote the Chern-Simons couplings) and 

■^*In the case oi k — 1, 2, the supersymmetry should get enhanced to d=3, 3^—8. 
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Table of Notation I 


Symbol 


Definition 


Symbol 


Definition 


Bulk S 


Dacetime 


9ab 

Gd+i 
h 

F{br) 

i 


Bulk metric with 
components u^, QJ^, 
Bulk Newton constant 
Inverse Horizon radius 

1 rf(br) 

2rDf(brD) 


Ad+i 
f{br) 
F{br) 


Bulk co-metric with 
components u'', qj^^QJ^ 
-ld{d-l) Bulk C.C. 
1 - {bry^ 

{F{br)-F{brD)) 


Boundary at r = oo 


V 
b 

c 

V 

f ^ 


Boundary metric [r~'^ds'^_^_-^]r=oc 

Fluid velocity 

Fluid pressure (D.17) 

Bulk viscosity (D.17) 
Christoffel covariant derivative 

— A p A 

Shear strain rate (A. 12) 
Acceleration field (A. 14) 


9^" 
p 

e 
s 
V 

T ^ 

e 


Boundary co-metric 

9iJiu + u^Uy 

Fluid energy density (D.17) 
Fluid entropy density 
^ Shear Viscosity (D.17) 
Christoffel symbols 
Weyl-Connection (A. 16) 
Fluid vorticity (A. 13) 
expansion rate (A. 15) 




Speed of sound-squared ^ = ^ 


Dirichlet Hypersurface at r = (Dressed dynamics) 


P 

a 

c 

V 

f ^ 


Dirichlet metric [i^~'^ds'^^-^]r=ro 

Fluid velocity 

Fluid pressure (A. 25) 

(1 - (brDy)-'/' = f-"\brD) 

Bulk viscosity (D.17) 

Christoffel covariant derivative 

— A F A 

Shear strain rate 
Acceleration field 


9^" 
P 

e 
s 

V 

t ^ 
9 


Dirichlet co-metric 

Fluid energy density (A. 25) 
Fluid entropy density 
£ Shear Viscosity (D.17) 
Christoffel symbols 
Weyl-Connection 
Fluid vorticity 
expansion rate 


^ D,snd 


radius where dressed sound becomes superluminal w.r.t g i.e., where Csnd = 1 


snd 


Speed of sound-squared ^ = ^ 1 + |(q;^ — 1) 
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Table of Notation II 



Symbol Definition 



Symbol Definition 



Incompressible Limit ala BMW(Reference [13]) (Add hats for fluid at r = td) 



-1 



W 

KM, 

V(o) 



BMW expn. parameter with p^,h^^dt ^ ^ ^ and v^, Vj ~ K 



-1 



(0) 



6o,ao, 



Po 



Spatial metric backgnd. 

Covariant derivative using glj 
Fluid velocity ^ ~ ^^"^ 
Backgnd. values 
Backgnd. pressure 
Pressure/mass density 



2^ ^ ~ 

^snd To 



PO 



ij 

K 

po 



Spatial co-metric backgnd. 

Vf ^fc* - Vfk* 
Mass density + Po 
Variation-typically ~ 
Backgnd. energy density 
Kinematic viscosity — 

^ _ rfjhqr) 



Incompressible Limit ala BKLS (References [15, 16, 19]) at r = r^j (with near Horizon) 



X 

6,, a 
p. 

P* 



BKLS expn. parameter with p^, dt ^ k 

Further 6,, d,, Csnd ~ k and h^, 



and , Vj ~ x 



-3 



(•) 



Spatial metric backgnd. 
Covariant derivative using g)' 
Fluid velocity |^ ~ 
Backgnd. values ~ x 
Backgnd. pressure ~ x^(^^^) 
Pressure/mass density 



snd 



* ~ X 



p. 



PD 



Spatial co-metric backgnd. 
Mass density 
e,+p, ^ x-^'^"^) 
Variation- ~ x 



-3 



Backgnd. energy density ~ x^'' 

Kinematic viscosity ^ 

p co-ordinate of Dirichlet surface 



P 



Near-horizon radial co-ordinate defined via r 



1 

X b. 



1+ 



a 't Hooft coupling A = Nc/k is conjectured to be dual to M-theory on AdS4xS|j._2/Zfc 
with G4 = Ar-2y9A78 = 3A;-i/2(2Arj-3/2. 

3. A d=6, K=(2,0) superconformal theory on M5-Branes is conjectured to be dual to 
M-theory on AdSyxS^^^/s with G7 = 37rV(16Ar3). 

The general bulk metric which is Weyl-covariant is given by 

ds"^ = QjiBdx^dx^ = -2 Uf.dx'' {dr + r ^^dx") + e^^dx^'dx" 
which is invariant under the boundary Weyl transformations 

where = 0(x) is an arbitrary function at the boundary. Without loss of generality assume 

0. Further we have 



(A.l) 
(A.2) 



that 0^1/ is transverse to u^, i.e. , ^^yU" 



Xl^[X) 



23. 



QJi,(r, x) 



(A.3) 
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We will raise/lower /contract the unhatted greek indices using the boundary metric 



g^ydx^dx" = [r '^ds^^ 



and the velocity field is a unit time-like vector of this metric 

u'^u^ = g^^n^Uu = -1 



(A.4) 



(A.5) 



A Weyl-covariant basis for the bulk cotangent bundle is given by {dr + r^^dx^ ,dx^} 
with a corresponding dual basis {dr , 9^ — r QJ^^r}. In this dual basis, the co-metric (or the 
inverse metric) is given by 

S^^Sa ® Ob 

= 2 uqS^ - r ^^dr] ®s dr + ^^'^ [d^ - r (g) [d, - r 2J,a,] 

where ^^'^'^ is the unique transverse tensor that satisfies 

*P^X = and ^^'^(5xu = S^ + u''Uy (A. 7) 



(A.6) 



The unit normal vector of a hypersurface r = is given by 
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HAdx^ 



dr 



dr 



_ (^"^ Q3„ - 2 r QJ„) dr + (u^ - QJ«) 9^, 

From these expressions, it follows that the normalized induced metric and co-metric on 
a hypersurface r = r^) is given by 



rv ^ 



(Ai 



The hatted greek indices are raised/lowered/contracted using this hatted metric/co-metric. 



■^^The corresponding problem for general function r = can be reduced to this problem by a suitable 

choice of Weyl frame. 
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A. 3 Notation at the Boundary 



The metric on tlie d dimensional boundary is denoted by gfj,^, which is a representative the 
class of metrics on the conformal boundary of the bulk spacetime. 



9f,u = {r "^ds"^}^^ 



(A.9) 



The inverse of this metric (we will call this the co-metric - since it is the metric on the 
cotangent bundle) is denoted by g'^'^. We denote with V the corresponding Christoffel con- 
nection/covariant derivative. Our conventions for Christoffel symbols and the curvature 
tensors are fixed by the relations 



V^V = d^V + r^A^V^-" and [V^, V.]V^^ = -R^uaV. 



(A.IO) 



On this spacetime lives a conformal fluid with velocity field (with u^u^ = —1) , 
pressure p, energy density e = {d — l)p, and shear viscosity rj. We introduce the projector 
Pfiu = dfiu + u^Uij which projects onto the space transverse to u^. 

The gradients of the velocity field are decomposed as: 

e 



P 



(All) 



d-l 



where we have introduced 

• the shear strain rate: 

• the vorticity: 

• the acceleration field: 

• the expansion rate: 



pa tdP 



P. 



I^i^ pa, 



d-l 



00. 



6 = V„m" 



(A 12) 

(A13) 
(A. 14) 
(A15) 



The hydrodynamic Weyl-connection (see [33] where it was introduced for more details) is 
defined to be 



-Ui, 



^ d-l 

The bulk metric-dual for hydrodynamics is given by (see [12]) 



(A.16) 



u 



&^u = P^.u + 2hF{hr)a^, + ... 



(6r)-^) u, + ... 
-- P^'' -2bF{br) a^"" + 



(A 17) 



where 



F{br) 



a 



{F{hr) - Fibro)) 



1 

(y. J ^j. 



bro yd-1 



ylyd- - 1 



-dy. 



(A18) 
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The energy-momentum tensor of a general relativistic fluid (till first order in the gradient 
expansion) is given as 

T^"' = {e + p) + p g^" - 2r] a^'' - C 9 P"" + . . . (A.19) 

which can be computed from the bulk data (of the metric dual to hydrodynamics) via 

T^, = ( "1 [-2 K^, + 2K g^, - 2 {d - 1) g^, + . . .]] (A.20) 



where Kfj^^y is the extrinsic curvature of the constant r hypersurface and K = g'^" K^i, is 
its trace. This computation gives 

^ 11 11,, , 

and the Bekenstein-Hawking argument in the bulk^° gives the entropy density and the tem- 
perature of this fluid as 

^ ^ and T = (A.22) 



4Gd+i 6^-1' A-Kh 
where we have found it convenient to introduce a variable h = (AGd+is)~'^ . 

A. 4 Notation on the Dirichlet Hypersurface 

We choose a hypersurface r = in the bulk to impose Dirichlet boundary condition. We 
will work in a boundary Weyl-frame where td is independent of x. This Weyl-frame change 
is consistent with the gradient expansion in the boundary hydrodynamics provided the initial 
surface r = p{x) had a slowly varying p{x). We will introduce a parameter 

a = ^ (A.24) 

which parameterizes how far the Dirichlet surface is to the boundary/ how close it is to the 
horizon. We have d(ro = oo) = 1 and d(r£) = 1/b) = oo. 

After this we proceed as we did in the r = oo case in the previous subsection. All 
the same definitions can be repeated - we will just distinguish the objects in the Dirichlet 
hypersurface by a hat - so we have u^, g'^", p and so on. Unless specified, all the hatted 
tensors are raised and lowered by the hatted metric/co-metric. 



^°Note that the above forms can be deduced upto a normahzation constant from just the thermodynamic 
relations 

ds de dT dp ^ ^ , . , , 

— = , — = — ^ and e+p = Ts A.23) 

s £+p T e+p ^ ' 

So the Bekenstein-Hawking argument is needed just to fix the normahzation constant. 
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The energy momentum tensor is calculated using the same expression as in equation 
(A. 20) except that we evaluate it at r = now. We get 



(A.25) 



where 



{d — 1) a 



2a 



P = 



V = 



Sit Gd+i ce + 1' ¥ a + 1 



a 



2a 



Sir Gd+i a + l ¥ a + l 
1 1 



d , 



P 



l&TiGd+i ¥- 



- = 7] and C = = C 



(A.26) 



B Dictionary for the Dirichlet Problem 

In this subsection, we collect the formulae which translate between the boundary data and 
the Dirichlet data to serve as a ready reference for the reader. The (normalized) metric and 
the co-metric on the Dirichlet hypersurface are given by 



1 - — ) UnUy + 2hF{hrD) (JaU ^ [UaAy + AnUy] + . . . 

a^ I To 



1 2^ \ 1 n , X 

+ I 1 - — + _ 1 UiiUv - — [u^,a^ + a^Uy\ + 2bF{brD) a^^ + ... 



g^" = Pf"" - 2bF{brn) cy^"" - 6? [-'^ 

2a^B 



a~ IM'* - Tj^a^^ \u'' - Tn^a 



9 



1-a^ - 



1 + 



a 



2a^9 



{d-l)rD\ 



, , u^'u" - 2bF(brD) a"^ + — lu'^a'' + a^'u''] 
{d-l)rDj Td 



and the correctly normalized velocities at the hypersurface are 

a^9 



Ufj, a 

Un = ~ I Aa 

a td 



rD{d-l)) a To 



Uf, a 



a^e 



rnid-l) 



+ 



(B.l) 



fB.2l 



From these it follows that the transverse projectors are related by 



P^.u = P^.u + 2bF{brD)a 



flU 



a 



a 



pi^u ^ p^u _ 2bF(brD)a'"' + — Ka'' + a^u" 



(B.3) 
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Given a relation of the form f/^ = V^, we can always write 



(B.4) 



which defines the tensor Tfj^^P = f ^jj,'' — T^^^. We evaluate this tensor in the Appendix C to 
get 



e 



d 



+ 



in particular 



r ''v 



{a' - 11 



From this it follows that 



(B.5) 



(B.6) 



a 



oj^y, 6 = a6, 



d 



a^, = + ^ {a^ - 1) a,. 



(B.7) 



We refer the reader to the previous subsection for the dictionary involving the energy- 
momentum tensor. 

Now we are ready to present the inverse relations. We will first invert the equations 
above to get 



a 



111/ ^ ^ riu ) ^riu 

a 



and then the velocities 



5 Au A-i/ 



a 



(d^-i); 



[l+|(d2-l 



(B.8) 



u,. 



1 + 



a9 



a 



rD{d-l) 



1 - 



a9 



a 



TD {d - 1: 



+ . . . 



(B.9) 



followed by the projectors 



2b 



Pfiu = Ppv- — F{brD) ^f,u 
a 



pp ^ 



a 



Td [I + - 1 



(B.IO) 
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The dictionary for the (normahzed) metric and the co-metric are 



rD(c?-l) 
1 29 



aroid — 1) 



rz) [1 + 1(^2-1)] a 
u^'u" + —F{brD)a^"' ^ ^— ^ 



(B.ll) 



Finally, we can write the tensor in terms of hatted variables as 



V = (1 - 



6 (i^2--^ da'^ 



+ {a' - 1) 



- 1 



In particular 



f^.% = -(l- 



(B.12) 



(B.13) 



Now wc present the Bulk metric/co-metric as a function of the Dirichlet data. The Bulk 
metric is given by 



SAsdx^dx^ = -2u^dx'' (dr + r A^ + ^{1- {br)-'^)u^ dx"^ 



where 



+ [P^^u + 2bF{br) a^,J\ dx^dx"" + . . . 



= -2 





. 1 


dx^ [dr + r 


2 


OLUa A a 







{dr + r iA,, + ^-{l-{br)-'^)au^ dx"^ 



P^, + 2bF{br) 



dx^dx" + ... 



A,, — 



" 2 

a r 



[l + f(d2-l)] 
(1 - (6r)-'^) 



t = 1 _ ^ ^ (1 - (6r)-'^) _ ^ _ 
F(br) = - (F(br) - Fibro)) = - / ^dy . 



(B.14) 



(B.15) 
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The Bulk co-metric is given by 



+ ^-2 [pM^ _ 2hF{hr) a^"] ® 



r\l - {hr)-") - 



2r9 



a{d — 1) 



dr (8) dr 



+ 2 



1 - 



(9^ ®5 9^ 



- 2bF{br) a 



In terms of the components in the Weyl-covariant basis, we have 



a 



d-i 



u,. 



= P,.u + 2hF{hr) a^, = P^, + 2bF{br) a^. 



+ ^{l-{br)-')&u. 



a 



a I 



rD{d-l) 



qjM'^ = pi"' - 2bF{br) a"" = P"" 



2bF{br) a^"" 



(B.16) 



(B.17) 



C Lorentz-Covariant derivative of the induced metric 



Wc wish to find the tensor V y^y'^ which describes the difference between the covariant deriva- 
tives at the boundary and the Dirichlet surface, i.e.. Given a relation of the form t/^ = V^, 
we can always write 



(C.l) 



which defines the tensor F^j,''. 



For definiteness, in this subsection we will continue to raise/lower/contract using the 
boundary metric gy^y - this means in particular raising/lowering/contracting do not commute 
with the hatted covariant derivative V so we need to be a bit careful. 
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As usual, zero-torsion condition implies T^,^'' = Fj,^'' and metric compatibility with g^j^i, = 
9t,y + h^y gives 

(C.2) 



where 



Since all our expressions are exact upto second derivatives, it is enough to work with just 
the zero derivative piece in h^^. 



a 



d2-l 



(C.3) 



It follows that 



r Pu 



ia' - 11 



(C.4) 



We can now evaluate 



P 

d-l 



u^a„ I 1 + -{a — r 



-(a - l)Af,Uy 



Finally, we obtain 



41 



fii/ 



a 



(C.6) 



from which it follows that 



I + -{&'- I] 



(C.7) 



= yi^ - (1 - a^)-a^ 



"^^We have used 
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We invert the last relation to get 



Jii, 



l-f(l-d2) 



d-1 



-u,. 



This can be used to write is terms of hatted variables 



a 



-J\,, 



ae 



a 



au 



rD{d-l) ^ rz5[l + |(d2_i 



(C.9) 



Now, we want to write Tuv'' in hatted variables. We start with 



T,/ = {a' - 1) 



da{^^Uu) 



+ 



6^-1 



and use the Dirichlet dictionary to get 



r 



da^ 



+ {&' - I) 



[1 + 



(C.IO) 



(c.ii: 



D Non-relativistic scaling a la BMW 

In this appendix we review and correct the scaling limit of [13]. The major change in nota- 
tion we make is to replace e-Quw by a parameter so that the BMW limit is a large K 
asymptotics of the expressions below. 

For simplicity we consider a fluid without bulk viscosity (which includes conformal fluids) 
with an energy-momentum tensor T'^'^ given by 



Tf^u = p Qf^u + + p) Uf,Uy - 27] + . 



(D-1) 



where we assume that this fluid lives on a background spacetime with a metric g^u and for 
the moment have just written out the stress tensor to first order in gradients. 



D.l Spacetime split for the non-relativistic scaling 

We begin by decomposing this metric into an ambient part g^l^} and a forcing part h^i,-, the 
split being done so as to recover explicit forcing terms in the Navier-Stokes (in addition to 
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the pressure gradient term). One picks a suitable frame for the ambient metric, and writes 
the geometry as 



= gf} + Ku , gj^J dx^ dx" = -de + g^^ dx' dx^ , (D.2) 



where g^j^ are slowly varying functions of with h^i, being treated as a perturbation. The 
metric perturbations which force the fluid are taken to be 

/i^^ dx^ dx" = 2 ^^-^ k* dt dx' + K'^ (/i*^ d^ + h*j dx' dx^) (D.3) 

where ^^ is the book-keeping parameter that implements the BMW scahng (note ^^ = ^bmw)- 

We employ the notation that all the functions which have a * subscript or superscript 
(which we freely interchange to keep formulae clear) are of a specific functional form with 
anisotropic scaling of their spatial and temporal gradients. 

y,(i,xO : W^-^'^ ^ M , such that {^^y^i, x'), Vf ^.(i, x^} - {0(^-'), 0(^-1)} (D.4) 

The co-metric corresponding to the metric above is given as 

gi^'^d^ ®d,^ -dt ®dt + g'l^-^di ®dj + 2 i^-^ ki dt ®s di 

- [{hi, - k*ki) dt®dt+ (h^ + Kki) di ® dj] (D.5) 
+ 2 [{h;, - k*ki) ki - hi^k*] dt^sdi + o{^-') 

where we have freely raised and lowered the spatial indices with g^j\ 
The velocity field of the fluid is parameterized as 

u" = u' (1,^-1 -y:) (D.6) 

where the function is determined via the constraint g^vU^u^ — —1. This gives the full 
velocity field in a large expansion as 



I.* = 1 + — (hi + 2 k* vi + ^j? vi vt) + 0{^-') 

= vl + ^ (hi + 2 k* vi + 4^ vi vl) vl + 0(K- 
^t^-l-\ {-h*tt + ^j? vi vt) + 0(^-') 



(D.7) 



Ui = {v* + k*) + 



Ki^i + \ {Kt + 2 k* vi + ^f,) vi vi) {vl + A;*) 



+ 0(R 



66 



which can alternately be written as 



^-2 

~2 



+ — + 2 k] vi + gf^ vi vt) dt + ^ {hi + 2 k] vi + gf^ vl v^) vld, + 0{^-^) 



^dx^ = -dt + {v* + k*) dx' - ^ (-hi + gf^ vi vi) dt 



Kjvi + ^ ( + 2 k* vi + gf^ vi vt ) {v* + k 



dx' + 0(K-^) 



(D.8) 



Now, we are ready to calculate the velocity gradients - with some foresight, we will use 
the fact that the BMW hmit is also an incompressibility limit where vl is divergenceless (see 
below). With this in mind, we can write the velocity gradients as 



e = 0(K- 



dtv* + vivfv* - /;j dx^ + 0(H-^) 

dtvl + viVf^l- f:]d, + Oi^-') 



a^^dxi'dx" = dx'dx^ - 2^-3 vi V^u^v*. dx'dt + Q{^-^) 



(D.9) 



a^^^d, ®d,^ ^-^Vl^vi^di ® dj + 2^-'{v* + k*)V^^^vi^dt ®s dj + 0{^-^) 
uj^^dx^ A dx"" = V[f t;*] dx' A dx^ - 2'^-\i "^fv*^ dx' Adt + 0(^-^) 

where V^-^'' is the covariant derivative compatible with g^'^\x'') and is the 'gravitational 
force' acting on the fluid 



/; 



d,k: 



vfk*-v';>k* 



7(0)7* 



dihl - dtk* + q*jvi 



(D.IO) 



where we have introduced g*- = V,-°^/c* — V^^^'k* 



D.2 Navier-Stokes equations on a curved geometry 

Now, we turn to the scaling of the thermodynamic variables. We define the mass density po, 
the pressure per mass density and the kinematic viscosity vq by 

Po = eo + Po , P = Po + ^~Vo P* and r/o = Po ^^o (D-H) 

where as before the subscript indicates the background value. All other thermodynamic 
variables have similar scalings, for example, 

e = 80 + Vo and 6 = 60 + (D. 12) 
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The BMW limit is taken as to be the scahng as K — )■ oo and in this hmit the conservation 
equation V^T'^'^ = reduces to 



and then a non-relativistic forced Navier-Stokes equation: 



0(K" 



dt + viV 



(0) 



(0) 



(D.13) 



(rf- 1) {d-2) 



where the kinematic viscosity z/q is given as 



1^0 



Po d 



(D.14) 



(D.15) 



Before proceeding we should explain the origin of this equation since it differs from 
that presented in [13]. On the l.h.s of (D.14) we see a familiar term corresponding to the 
convective derivative of the non-relativistic velocity. The usual Laplacian term is modified 
due to the background curvature into the second derivative piece multiplying z/q. Its origin 
can be traced back to the term —2iria^i, in the relativistic stress tensor (D.l). On the r.h.s. 
of (D.14) we have collected all the forcing terms: there is the familiar pressure gradient term 
along with two other terms that arise from curvature. /* is the forcing term that arises from 
the fluctuating part of the metric as is clear from (D.IO); this term has been accounted for in 
[13]. However, we also should see a forcing of the fluid from the 'background' curvature: the 
spatial part of the metric (]\iv is a curved spatial metric g^^^ and its effect on the fluid turns 
out to be at the order 0(K^^), just the same as the other terms in the equation. However, 
its origins in the relativistic stress tensor are a bit more involved; it does not arise from 
any of the terms written down in (D.l) but rather from a second order gradient term in the 
relativistic stress tensor 2rjhC^aui3u"'u^ i.e., a coupling between the fluid and background 
curvature [12]. 

Note that we have here specialized to relativistic conformal fluids which have holographic 
duals, so that the transport coefficient multiplying the tensor structure C^aupu'^u^ is fixed 
to be 2rih. For a general fiuid we can have a new transport coefficient here k oc rib and 
the correspondingly we would replace the coefficient of i?^"^ in (D.14) with Uq — )■ Also, 
we should note that that while other tensor structures involving curvatures couplings are 
allowed for non-conformal fluids, these will necessarily have non- vanishing trace and as a 
result will only show up at sub- leading order in the BMW scaling limit. 

At the risk of being overly pedantic we reiterate the fact that if we wish to place consider 
the non-relativistic BMW scaling of a relativistic fluid, then we must necessarily work with 
higher order gradient terms in the relativistic fluid stress tensor. To obtain the correct 
non-relativistic equations up to the order where we encounter Navier-Stokes equations the 
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relevant part of the relativistic stress tensor is given to be 



(D.16) 



which is a subset of the full second order stress tensor derived in [12] 



T^,y = pg^iu + {.e + p) u^u^ - 2r]a^^ 



V 



X 



(D.17) 



+ 27]b u^Vxa^^ + a^^axu- 



d-l 



It is a simple exercise to verify that none of the other terms involved in the second order 
stress tensor (D.17) contribute to the BMW scaled equations at 0(K~^). 

D.3 The bulk metric dual to a non-relativistic fluid on the boundary of AdS 

One can also construct the gravitational solutions dual to such fluids as described in [13]. 
To do so we simply need to apply the scaling of parameters described earlier as for e.g., 
in (D.6), (D.12) to the general fluid/gravity bulk metric dual to a relativistic fluid on the 
boundary of AdS. Such a metric correct to second order in the relativistic gradient expansion 
was originally derived in [12] generalizing the original result of [8]. It was believed that this 
in general would suffice to find the gravity dual for a non-relativistic fluid that satisfles the 
incompressible Navier-Stokes equations derived above (D.13), (D.14). In fact the original 
computation presented in [13] argued that it would actually suffice to consider the relativistic 
metric accurate to first order in gradients. Unfortunately, this turns out to be incorrect and 
one needs a subset of second order gradient terms along with one particular third order 
gradient term to solve Einstein's equations to order 0(K^^). Note that this is necessary 
because it is at 0(K~^) that we encounter the dynamical content of the boundary fluid 
equations, viz., the Navier-Stokes equation (D.14). 

The new ingredient in our analysis is that we need to worry about a particular third order 
term proportional to 'D^'Jl in the fluid/gravity correspondence. The term in question turns 
out to be computable using the original algorithm outlined for constructing bulk metrics 
dual to boundary fluids in [8] and thankfully involves a decoupled tensor structure that can 
be sourced independently. Including this term, we flnd that for the non-relativistic fluid 
on the Dirichlet surface it suffices to consider the following truncation of the relativistic 
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fluid/gravity metric to third order in boundary gradient expansion: 
ds^ = —2u^dx^ {dr + r Aydx") 



+ 



3r (rf- l)(ci-2) 



- 2 (for) 



H,[br)C,^,pU u + (rf_i)(rf_2) 



(D.18) 



+ ... 

where the functions appearing above and their large r asymptotics are given to be: 

1 



/(for) 
F(for) 



dy 



L{br) = 



{brY 

yd-1 _ I 

y{y^ - 1) 

^ ^ 1 # 

for d{brY ^ {d + l){brY+^ ^ {brf'^ ^ ' 

^'-'du dy y 



br 



1 

+ 



1 



(ci+l)(for) 2(rf + 2)(for)2 {d + l){2d + l){brY+^ 
1 , # 



(rf + l)(2rf + 4)(for)'^+2 (for)2<i+i 



+ 



(D.19) 



iJi(for) = 



00 yd^2 _ I 

y{y^ - 1) 

1 1 



1 



2(for)2 d{brY {d + 2){brY+^ {brY^ 



+ ... 



iV(for) = / e'^-irfe / dy 



br 



y\y^ - 1) 
1 

+ 



1 



(ci+l)(for) 3(rf + 3)(for)3 ' (rf + l)(2rf + l)(for)'^+i 
1 # 



+ 



(d + 3)(2d + 3)(for)'^+3 (for)2'i+i 



+ . . . 



There are various new curvature tensors and derivatives introduced above. 2) denotes 
the Weyl covariant derivative introduced in [33] which was used to present the bulk metric 
dual relevant for fluid/gravity correspondence in the case of curved boundaries in [12]. "R^y 
is likewise a Weyl covariant Ricci tensor and S^,y is a Weyl- covariant Schouten tensor 



1 



d~2 
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Qfiu "R 



(D.20) 



For further details of these objects and the complete form for the second order fluid/gravity 
metric accurate to second order in gradients we refer the reader to [12]. 

Once we have the relativistic metric at hand it is a simple matter to employ the scalings 
outlined earlier. We find that the bulk metric dual to an incompressible Navier-Stokes fluid 
living on a spatially curved geometry at the boundary of AdS takes the form 

ds"^ = dsl + ^-^ds\ + ^-"^dsl + ^-^dsl + 0(K-^) 
with 



dsl 



2 dt dr + 



-fodt^ + gifdx'dx' 



dsl = -2 « + k*) dx' dr + 2r^ [k* - (1 - /o) {v* + k*)] dx'dt 



dsl 



2"'tt ^ 2^jk 



dt dr + [/i*^ dt^ + h*j dx' dx^] 

+ r'{l- /o) [(^-h;, + gf^ vi vt + V*d) dt' 

+ {v* + k*) {v* + k*) dx'dx'] +2r' feoi^oVjf w*) dx'dx^ 



{d-l){d-2) 



dt' 



2H. S 



(0) 



i?(0),(0) 



2(c/- 1) {d-2) 



doc dx 



{D.21] 



dsl 



-2 

+ 2r 



dtv* + vivf\* - f* dx'dt - Ar%FQvlV';','v*. dx'dt 



[vl + k* 



(0). 



dx^dr 



2 (1 - /o) irl^vi + ( k* vi + gf;} vi v', + pj) {v* + k*) 

T-r2 J. A o o(0) 



-2 



+ 



dx'^dt 



j^^%)v: dx^ dt - 2 ^ vi dx^ dt-^ V^^g*^ dx^ dt 



+2 



{d-l){d-2) ' ' 
bo No d-3 



k*)dx'dt + AHo 
Vf^R^^Ux'dt 



(0) 



i?(0)^(0) 



2(d- 1) {d-2) 



vi dx" dt 



where we have highlighted the terms that were missed in the previous analysis for quick 
comparison. Note that when the background spatial metric gf^ is Ricci fiat, many of the 
terms vanish except for two terms in the third order metric (which are proportional to 
V(Q)f* and V(Q)g*j respectively). Finally we should note that S-^"* is used to denote the 

spatial components of the Schouten tensor of the full background metric g'^^}] in particular, 
it should not be confused with the Schouten tensor of the spatial metric gf^\ The functions 
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that enter into the metric above are: 



/o = 1 - {hry^ , = — and Fq = / -— —dy 

bo Jbor y{y - 1) 



y-i 



poo d-2 _ 1 

Ho ^ {h,rf / , ^ dy 



Jbor y^y - 1) 

j-oo POO 2 _ 1 



'bor 

Finally, let us note that the Navier-Stokes equations themselves have an interesting scal- 
ing symmetry. Given any solution to (D.13) and (D.14) with gl,j^ = 6ij we can consider 
replacing 

P*^e'p«, vl^evl, fl^f^fl (D.23) 

where again the functions entering the dynamics with subscript *e have spatial gradients 
(9j ~ e and temporal gradients 9t ~ e^. This fact makes it possible to compound the Navier- 
Stokes scaling which effectively allows one to replace K — )■ K'" for some w > 1. Essentially 
the incompressible Navier-Stokes system of equations is a fixed point set of this scaling 
symmetry, a fact that we have made use of in §6. 



E Bulk dual of the non-relativistic Dirichlet fluid 

In this appendix we present without derivation the results for the non-relativistic scaling 
limit of the Dirichlet problem, generalizing the result quoted in §5.2. Physically the only 
new content is that we allow the metric on the Dirichlet surface to be endowed with 
an arbitrarily slowly varying spatial metric. Thus in contrast to §5.2 we are relaxing the 
constraint on gf^ introduced in (5.2) being Ricci flat. 

To indicate the differences note that the presence of a non- Ricci fiat spatial metric gf^ im- 
plies that the non-relativistic metric gets contributions from various tensor structures which 
appear at the second order in the gradient expansion of the fluid/gravity correspondence. 
In terms of the metric written down in [12] some of these are straightforward to see - any 
term involving curvature tensors of boundary data (and thus via the Dirichlet constitutive 
relation the hypersurface S^i curvatures) will contribute at this order. However, we also get 
contribution from spatial gradients of the hypersurface curvature, i.e., in the BMW scaling 
limit encounter terms of the form To guide the reader towards a derivation, we 

quote simply the relativistic tensor structures which are relevant and their scaling behavior 
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under the Dirichlet BMW scaling. 



dx^ dx" = -de + gf^ dx' dx^ + 2 k* dt dx' + K"" ( /i*^ di' + h*. dx' dx^ 



u 



= -dt + [v* + rfx^ - ^ (-/i,*, + vi vi) dt 



+ \ {kt + 2 k* vi + gf^ vi v^) 



v: + k: 



d^dx^ = 



dtv* + t)^ Vf {)* - /; dx' + 0(K- 



a^^dx^'dx" = Vjf {)*) rfx^dx^' - 1)^' vjf {)*) rfx^dt + 0{k 
as before along with new tensor structures: 



S^xu^dx" = 



Jl^xu^dx" = K 



2(rf- l){d-2 



P^^KV -L- -\7^ n* 



^ij ^* + r,/ J ON ^coiyj 



2{d - 2) (0)^'^; 



dxi' + 0(K-^) 



-%~^v:dx' + o(s;^-^) 
MO) _ m 

2(rf-l)(d-2) 



-(0) 



dtij doc 



2K 



-3 



(0) 



a(o) „ 

2{d - l){d - 2f'' 



vidx^dt 



where 



(E.2) 



vf^k* - vfk* 



AKt - dtk* + q*jVi 



(E.3) 



(0) 



and we use S^^ to denote the spatial (i.e. ij-) components of the Schouten tensor for gi^J 
keep the expressions somewhat compact. 



E.l Boundary data for non-relativistic fluids on Tjd 

Given the scalings to obtain the non-relativistic fluid on S/), it is possible to identify the 
relevant terms of the third order fluid/gravity metric that we need to retain to solve the 
Dirichlet problem. Per se the set of terms we need in the bulk metric is still given by 
(D.18), though as in the main text we want to use this information to solve for the Dirichlet 
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constitutive relations. We first quote the results for the hypersurface stress tensor which 
defines the hypersurface velocity before indicating the answers for the boundary velocity 
field and metric in terms of the hypersurface data. 

We can parameterize the hypersurface stress tensor as in the main text; to obtain the 
correct non-relativistic equations on E/) we need to retain some second order gradient terms 
involving hypersurface curvature tensors (analogously to the situation at the boundary as 
described in §D.2). The relevant piece of the relativistic hypersurface stress tensor turns out 
to be: 



Tf^t, = {i + p) Uf,Uy + p g^u - 2fi a + kc C ^c,upu"'uP + . . . 



(E.4) 



where 

e 

i + p 
V 



d-1 



a 



da 



a 



2rl {d-l){d-2) 



1 
1 



2rl {d-l){d-2) 



a 



8nGd+i¥a + lrl{d-l){d-2) (g.S) 



1 - 



d-2 



a + l (bro) 



Conservation of this stress tensor together with the scaling forms introduced in (E.l), (E.2) 
leads to the incompressible Navier-Stokes equations on E^i. 

The boundary velocity field and metric can be expressed in terms of the Dirichlet data 
as before. Before we write out the exressions in their gory detail, let us introduce some new 
parameters k^, kjy which depend on the location of E^i as 



K-N = ^ {d - 3) 



d-2 



2^ [1 + f(d2- 1 



^=bri 



i=hrB 



(E.6) 



These are in turn limiting values of certain functions which appear in various guises in the 
result for the Dirichlet constitutive relations and the bulk metric and are collected once and 
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for all below. 
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(E.7) 



2b^[l + l{a^-l) 

where the functions entering the above expressions and their asymptotics have been previ- 
ously been collected together in (D.19). 



The final result of this exercise leads to: 
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Further the Dirichlet constitutive relation for the boundary metric is 
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E.2 The bulk dual for arbitrarily spatially curved metric on S/j 



The final result for the bulk metric dual to the non-relativist ic fiuid living on the Dirichlet 
hypersurface Ed is simply obtained by plugging in the scaling form (E.l), (E.2) into (D.18), 
having eliminated the boundary data in favor of the hypersurface data using the Dirichlet 
One obtains: 

ds^ = dsl + ^'^dsl + ^~^dsl + ^-^dsj + 0(K-^) 

with 
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